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FOREWORD 


This  report  was  prepared  by  the  Douglas  Aircraft  Company,  Long  Beach, 
California,  for  the  Aerospace  Dynamics  Branch,  Vehicle  Dynamics  Division,  Air 
Force  Flight  Dynamics  Laboratory,  Wright-Patteraon  Air  Force  Base,  Ohio,  under 
Contract  F33615-70-C-1167.  This  research  was  conducted  under  Project  1370, 
"Dynamic  Problems  in  Military  Flight  Vehicles,"  and  Task  137003,  "Prevention 
of  Dynamic  Aeroelastic  Instabilities  in  Advanced  Military  Aircraft."  Mr.  S. 

J.  Pollock  of  the  Aerospace  Dynamics  Branch  was  Task  Engineer. 

This  report  consists  of  two  parts  with  two  volumes  for  each  part.  This 
volume,  Volume  II  of  Part  I  is  the  Computer  Program  H7WC.  Volume  I  of  Part  I 
consists  of  the  method  of  direct  application  of  nonplanar  lifting  surface 
elements.  Volume  I  of  Part  II,  contains  a  method  which  uses  an  image  system 
and  an  axial  singularity  system  to  account  for  the  effects  of  the  bodies  and 
Volume  II  of  Part  II  is  the  Computer  Program  N5KA. 

The  work  reported  herein  was  conducted  during  the  period  of  December 
1969  to  September  1970. 

The  Principal  Investigator  was  Joseph  P.  Giesing.  Mrs.  Terez  P.  Kalman 
was  responsibxe  for  implementing  the  method  on  the  computer.  Donald  H.  Larson 
aided  in  this  implementation  and  acted  as  consultant  for  computer  problems. 

Dr.  William  P.  Rodden,  a  McDonnell  Douglas  Company  Consultant,  contributed 
many  valuable  ideas  to  the  project.  Others  have  made  significant  contribu¬ 
tions  to  this  project  including  Messrs.  D.  S.  Warren  and  W.  E.  Henry. 

The  report  was  submitted  by  the  authors  in  November  1970  for  publication 
as  an  AFFDL  Technical  Report. 
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ABSTRACT 


Two  methods  of  accounting  for  body-liftir  interference  in 

unsteady  flow  are  considered.  The  first  method  is  described  it.  Fart  I 
of  this  report,  while  the  second  will  be  described  in  Part  II  to  follow. 

The  first  method  is  a  direct  application  of  nonplanar  lifting  surface 
elements  to  both  the  lifting  surfaces  and  the  body  surfaces.  The  body  is 
treated  as  an  annular  wing.  This  type  of  idealization  must  be  used  with 
an  axial  doublet  introduced  to  account  for  body  incidence  effects.  The 
undesirable  effects  of  the  annular  wing  representation  are  then  reduced. 

The  second  approach,  to  be  described  in  Part  II,  uses  an  image  system 
and  an  axial  singularity  system  to  account  for  the  effects  of  the  bodies. 

This  volume  contains  the  computer  program  and  contains  the  FORTRAN 
listing.  Volume  I  contains  the  development  of  the  theory,  correlation  of 
theory  with  experimental  data,  and  the  parametric  study. 
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1.0  INTRODUCTION 


Computer  program  H7WC  is  the  result  of  modifying  and  extending 
program  H7WB.  Piogram  H7WB  generates  aerodynamic  influence  coefficients 
and  gust  loads  for  nonplanar  surfaces.  The  modifications  to  this  program 
include: 

(1)  An  extension  to  accept  polynomial  mode  input.  (Subroutine 
GENQ) 

(2)  The  addition  of  a  capability  to  determine  generalized  forces 
(Subroutine  GENF) 

(3)  The  addition  of  a  capability  to  account  for  body  angle  of  attack 
and  camber  effects  using  an  axial  pressure  doublet  system. 
(Subroutines  AUGW,  TKER) 

(4)  A  refined  version  of  the  Doublet -Lattice  Method  to  properly 
handle  problems  involving  small  vertical  gaps  between  wing 
and  tail  surfaces.  (Subroutine  IDF2) 

(5)  The  addition  of  the  Quasi-Inverse  for  sc  Ang  large  numbers 
of  modes  or  submodes,  and  for  solving  future  modeo  with  a 
minimum  of  computing  effort.  (Subroutines  QUAS,  FTITSOL, 

CXSS) 

Program  H7WC  can  handle  almost  any  type  of  configuration  if  it  is 
idealized  properly.  The  axial  pressure  doublets  give  the  proper  slender  body 
flow  field  for  an  axially  symmetric  body  of  any  shape.  The  body/lifting-surface 
interference  is  taken  care  of  using  lifting  surface  panels  on  the  body  surface 
The  interference  elements  must  lie  on  an  idealized  body  of  constant  cross 
sectional  shape.  The  rules  that  apply  to  the  placement  of  a  grid  of  boxes 
on  the  aircraft  lifting  surfaces  also  apply  to  the  interference  surfaces  on 
the  bodies.  These  rules  are: 

(1)  Boxes  (trapezoidal  elements)  are  arranged  in  strips  parallel 
to  the  free-stream.  The  aspect  ratio  of  such  boxes  should  not 
be  large.  For  the  unsteady  case,  an  aspect  ratio  of  order  unity 
or  less  is  preferred. 
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(2)  Surface  intersections,  surface  edges,  fold  lines,  and  hinge  lines 
should  lie  on  box  boundaries. 

(3)  For  wing-tail  type  configurations  the  strip  boundaries  on  the  tail 
must  be  aligned  with  those  on  the  wing.  Strips  should  be  concen¬ 
trated  near  the  wing  tips,  near  control  surface  edges  and  in  any 
region  where  the  span  loading  changes  rapidly. 

(4)  The  box  length  must  be  small  compared  to  the  basic  wave  length. 
— •  <  0.  04 

Boxes  should  also  be  concentrated  near  control  surface  hinge-lines 
and  in  all  regions  where  the  upwash  boundary  conditions  are 
discontinuous. 

(5)  The  body  cross-sectional  shape  on  which  interference  lifting 
surfaces  are  placed  must  be  constant  in  the  x-direction  even 

for  wing-tail  problems.  (This  does  not  apply  to  the  slender  body 
solution  for  which  the  radius  distribution  is  used.  ) 

Sketch  1,1-1  gives  examples  of  these  5  rules. 

Because  of  the  numerical  techniques  involved,  the  two-dimensional 
radial  distance,  J  ( y-%)2  +  (z-4e)2  •  between  a  control  point,  y,  z.  (3/4-chord 
point  of  a  box)  and  a  box  edge,  ne>  4e>  may  never  approach  zero.  It  does 
not  matter  whether  the  control  point  is  on  the  same  or  on  a  different  sur¬ 
face  orwhether  it  is  up-or  downstream  from  the  box  in  question.  This  is 
easily  understood  for  a  control  point  in  the  wake  of  the  box  since  the  con¬ 
trol  point  may  not  approach  the  bound  leg  of  th e  horseshoe  vortex  However, 
it  is  also  true  upstream  due  to  the  approach  taken  for  the  unsteady  effects. 

(If  k  =  0  the  upstream  effect  at  o  (y-H)^  +  (Z_C)  -*  t)  vanishes). 

Wing-tody-tail  problems  must  be  hanuied  using  only  one  interference 
cross-section.  That  is, one  constant  cross-sectional  shape  must  be  used 
both  near  the  wing  and  near  the  tail.  If  this  rule  is  not  followed  the  wakes 
of  the  lifting  surface  element  lying  on  the  c  ros  s -section  near  the  wing  may 
impinge  on  the  tail  surface. 
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When  a  body  mode  is  input  it  must  be  input  both  to  the  axial  elements 
and  to  the  interference  body  panels.  It  must  be  remembered, however, that 
the  body  is  in  motion  either  in  the  z-or  y-directions  while  the  lifting  surfaces 
on  the  body  are  in  motion  normal  to  themselves.  Thus  the  mode  input  for  a 
particular  body  interference  panel  must  be  multiplied  by  cosVfor  (+)  z-motion 
and  -sinY  for  (+)  y-motion  where  Y  is  the  dihedral  angle  of  the  panel. 

The  following  list  gives  the  program  limits. 

1.  The  maximum  number  of  unknowns,  i.  e.  ,  the  aerodynamic  boxes 
on  the  lifting  surfaces  plus  the  number  of  body  elements,  is  400. 

2.  The  maximum  number  of  aerodynamic  degrees  of  freedom,  in 
case  of  AIC  calculations,  is  130. 

3.  The  maximum  number  of  degrees  of  freedom  per  strip  in  case 
of  AIC  calculations  is  restricted  to  4.* 

4.  The  maximum  number  of  modal  right-hand  sides  is  40  per  case. 

5.  The  maximum  number  of  aerodynamic  boxes  per  spanwise 

division  (strip)  is  30. 

6.  The  maximum  number  of  spanwise  divisions  (strips)  is  50. 

7.  The  maximum  number  of  bodies  is  20. 

Program  H7WC  was  used  successfully  on  the  following  computers:  GE635. 
IBM  7094  and  IBM  360/65.  Core  requirement  for  the  program  on  the  various 
computers  is  as  follows: 

CE635  IBM  7094  IBM  3bO/65  CPC  6600 

28K  words  32K  words  12f>K  bytes  I00K  bytes 

If  greater  core  capacity  is  available  it  is  desirable  to  use  it  up  to  2b0K  bytes. 
The  CPC  version  is  slightly  modified  in  dimensions. 

*  This  restriction  comes  timn  operational  requirements  on  the  IBM  7Q'M 
The  number  of  degrees  of  freedom  can  easily  be  increased  to  the  7 
discussed  in  Appendix  P  in  Part  I,  Vol  I  t>y  increasing  ti  e  dimension  4  to  7 
in  subroutines  Alt',  AERO,  and  FINAL. 
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2.0  COMPUTER  PROGRAM  H7WC 

2.1  Input  Procedure 
2.  1.1  Input  Sheets 

The  input  sheets  for  program  H7WC  are^hown  on  the  next  three 
pages.  The  first  three  cards  represent  general  ca^e  data  that  is  input  once 
per  case.  The  next  four*  cards  (4,  5,  6,  7)  represent  panel  data  that  is 
repeated  per  panel.  The  panel  data  must  be  inprt  in  a  certain  order. 

The  panel  data  for  all  "regular"  lifting  surface,  (wing,  tail,  etc.  )  are 
input  firs>..  The  panel  data  for  all  interfere>ce  lifting  surfaces  (used 
for  the  bodies)  is  input  last. 

If  a  body  can  oscillate  in  bon  th e  z-  and  y-directions 
(nacelle)  then  it  must  be  replaced  by  K0  bodies:  one  that  moves  in  the 
z -direction  and  one  that  moves  in  tN  y-direction.  The  subsequent 
three**  cards  (8,  9.  10)  are  body  d,ta  necessary  for  the  slender  body 
analysis.  This  data  is  repeated.jer  body.  Card  11  is  used  for  mis¬ 
cellaneous  modal  and  AIC  dataand  must  be  completed  in  euher  case 
Card  12  is  used  when  a  gust  s  desired.  Card  13  indicates  a  series  of 
cards  that  are  used  to  iden*iy  lifting  surface  boxes  with  strips.  Card  14 
indicates  one  or  more  car18  that  identify  the  degrees  of  freedom  to  be 
used  for  each  panel.  Ths  card  is  necessary  only  if  AIC's  are  desired. 

Cardsl5  a* d  16  are  used  to  identify  axes  for  the  AIC  degrees  of 
freedom.  Card  17  s  used  for  modal  control  data.  Card  18  is  used  to 
input  control  dataxor  Panel  modes,  while  card  19  represents  a  series  of 
cards,  repeater’as  necessary,  to  indicate  polynomial  mode  coefficients 
for  panels.  CirUs  4Q  a»d  21  are  similar  to  cards  18,  19  except  the  modes 
pertain  to  t’a  z  anc*  Y  modes  of  the  bodies. 


-here  may  be  more  than  four  cards  depending  on  the  number  necessary 
to  present  all  0  and  r  values 

i  i 

There  may  be  more  than  three  cards  depending  on  the  number  necessary 
to  present  F  and  RAD. 


5 


DO  NOT  PUNCH  RLANK  COLUMNS  _  NO  UNDE RPUNCHE S IN  SIGN  FIELDS. 


REPEAT  PER  PANEL 


VtCmtH _ _ _ PHONE  _  DATE 


OlffCCriONl  FO*  KEY  PUNCH  DO  NOT  PUNCH  BLANK  COLUMNS _ _  NO  UNDEftPUNCHES  IN  SIGN  FIELDS. 


2.  1.2 


Description  of  Input  Data 


Item  No.  Description 

1  The  header  card  is  used  to  describe  the  current  run  being  made. 

Any  alphanumeric  information  can  be  used,  placed  in  cc  1-60. 
The  first  12  card  columns  should  contain  the  tape  title  for  the 
AlC's,  in  the  case  when  AIC's  are  generated  and  are  to  be  saved 
on  master  tape  for  later  use. 


2  M,  FMACH  —  Mach  number,  usual  definition 

3  A,  ACAP  —  Reference  area;  usually  total  area  of  wing 

4  c ,  R E F C H D  —  Reference  chord,  usually  average  chord  of  wing 

5  s,  REFSPN  —  Reference  semispan. 

6  6,  NDELT  —  Symmetry  flag  ( 6  =  1  sym.  ,  6=  -1,  antisym 

6  =  0,  no  sym.  ). 

7  NP,  NOPAN  —  Total  number  of  panels  on  all  lifting  surfaces 

8  NB  —  Total  number  of  bodies. 

9  NRF  —  Number  of  reduced  frequencies  for  which  the  analysis 

has  to  be  performed  (repeated)  -  recommended  maximum 
is  6  per  case. 

10  NCORE  —  Size  of  total  solution,  i.  e.  ,  NCORE  N  x  M  where 

N  is  the  total  number  of  elements  or  unknowns 
(including  z  and  y  fuselage  elements)  and  M  is  the 
number  of  modes  or  the  order  of  the  Alt'  matrix  if 
obtai ned 


11  N !  —  If  1".  program  assumes  that  case  has  been  run  previously 

and  that  a  quasi -inve  rse  exists  on  a  tape  (identified  by 
NIP!)  and  only  new  modes  are  being  run  it  "0  ,  program 
runs  a  new  case. 
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Item  No. 


Description 


12  N2  —  If  "0", polynomial  modes  are  to  be  input.  If  "1"  AIC's 

are  to  be  calculated. 


13  N3  —  Data  flag.  If  "1",  matrix  is  printed;  if  "0",it  is  not  printed. 

14  N4  —  Detail  print  flag  for  slender  body  solution;  1  to  print,  6  otherwise. 

15  KR  —  Array  of  reduced  frequency  .  KR  =  . 

L  Uqo 

16,  17,  18  XI  X2  X3  —  Panel  edge  coordinates 

19,20,21  X4Y1Y2  XI,  Yl,  Z 1  (Inboard  leading  edge) 

22,23  Z1  Z2  X2,  Yl,  Z1  (Inboard  trailing  edge) 

X3,  Y2,  Z2  (Outboard  leading  edge) 

X4,  Y2,  Z2  (Outboard  trailing  edge) 

24  NS  —  Number  of  spanwise  divisions  for  panel  considered 

25  NC  —  Number  of  chordwise  divisions  for  panel  considered. 


SKETCH  2.  1.  2-1 
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Item  No.  Description 

26  COEF  —  Scale  factor  for  panel  deflection  modes.  If  no  scale 
factor  is  desired,  must  be  set  to  1.0. 

27  G  ~  Chordwise  divisions  in  fraction  of  chord.  Usually  varies 
from  0  at  leading  edge  to  1 .  0  at  the  trailing  edge. 

28  t,  —  Spanwise  divisions  in  fraction  of  panel  span  Usually 
0.  0  at  inboard  chord  and  1 .  0  at  outboard  chord. 

Repeat  all  panel  data  per  panel  (Items  16  through  28)  Panels  are  input  in  the 
following  sequence-  wing,  taibetc.,  panels,  then  body  interference  panels. 

29  —  Z-coordinate  of  body  axis. 

30  —  Y  -coordinate  of  body  axis. 

31  NF  —  Number  of  divisions  on  body. 

32  NZ  —  If  "1",  the  program  will  allow  for  doublets  in  the 

Z -direction  and  thus  will  be  able  to  account  for  Z-or  upwash. 

If  not,  set  NZ  0. 

33  NY  —  If  "1",  program  will  allow  for  Y -doublets  and  thus  Y-or 
sidewash . 

34  COEF  —  Scale  factor  for  body  deflection  mode 

35  MK1~  Each  body  has  associated  with  it  interference  lifting  sur¬ 
faces.  MK1  represents  the  first  box  on  these  surfaces 

36  MK2-  The  last  box  on  the  bodv  interference  lifting  surface 
panels.  Interference  lifting  surface  panels  are  input 
after  the  regular  lifting  surfaces  are  input. 

37  £  -  Endpoints  of  body  divisions  (starts  at  leading  edge  and 

ends  at  trailing  edg'-l.  F  x  coordinates  of  bodv  element  end  points. 

3  8  RAD  —  Radius  at  entl points 

Repeat  all  body  dsta  per  body  Uterus  2'1  through 


SKETCH  2.  1 . 2-2 


Desc  nption 

NSTRIP  —  Total  number  of  strips  on  all  lifting  surfaces  . 

NB1  -  Sequence  number  of  first  strip  Omit  if  no  AIC's 
are  desired . 


1  N B2  —  Sequence  number  of  last  strip  Omit  if  no  AH'  s 
are  desired. 

2  NPR1  —  Print  flag.  If  \",  print  all  AIC'  solutions.  Otherwise 

use  "0"  (usually  -  0). 

1  N  PR  2  —  Print  flag.  If  "  1 ",  print  AIC's.  Otherwise  use  0 

(usually  0) 

4  NPR1  —  Stability  derivati.e  flag  If  "I",  prim  static  stability 

derivatives  If  "  4" .  print  s  tat  it.  -  and  dynamic  stability  denva'.tes 

Otherwise  use  0  Note  that,  if  NPR  '  v  NRF  mus’  be  '  2  ami 

FR(1)  must  Vie  0  0  i'suai  setting  ;s  Ni'K  '»! 


OB  NO  —  AH'  tape  identify  at:  »n.  Ign.  re  if  no  AH"  s 
.1 QBNQ2  -  Net  used 

N(,rs;  -■  Oust  flag  i‘s<-  !  for  gust.  0  otherwise 


item  No. 


Description 


4C  JSP  —  Second  symmetry  plane  (the  plane  Z  =  0). 

-1  ground  effect  (antisymmetric) 

1  biplane  or  jet  (symmetric) 

0  no  symmetry 

49  NPC  —  Mode  selector  for  A1C  generation. 

Ignore  if  no  AIC's.  Use  "0"  if  alternative  #1 
is  used  i.  e.  plunging,  pitching,  control  sur¬ 
face  rotation  and  tab  rotation.  Use  "1"  if 
alternative  #2  is  ^sed.i.e.,  three  cambering 
modes  plus  control  surface  and  tab  rotation. 

50  NSV  —  Total  number  of  strips  on  all  vertical  panels 

which  lie  on  the  symmetry  plane  Y  =  0. 

Vertical  panels  lying  on  Y  -  0  must  be  input 
before  the  other  panels.  For  example  the 
vertical  portion  of  the  T-tail  must  be  input  first. 

51  NBV  —  Total  number  of  boxes  on  vertical  panels. 

52  NY  AW  -  Stability  derivative  override  flag  If  "1",  yaw 

coefficients  will  be  calculated,  if  "0",  pitch 
coefficients  will  be  calculated. 

If  NDELT  =  1 ,  NY  AW  =  0 

NDELT  =  -1,  NY  AW  =  1 
NDELT  =  0,  NY  AW  =  0  or  1 

53  GZRO  -  Gust  reference  plane  dihedral  in  degrees.  For  0 

dihedral  (GZRO  =  0)  the  gust  velocities  are 
vertical  (in  the  Z -direction) .  Omit  if  gust  is  not 
desired. 

54  XZRO  -  Gust  reference  point  Point  at  which  sinusoidal 

gust  velocity  is  unity.  Omit  if  gust  is  not  desired. 
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Item  No. 


Description 


55  VEL  —  Aircraft  /elocity  in.  / sec.  Omit  if  gust  is  not  desired  • 

56  WG  —  Gust  vertical  velocity  at  XZRO,  in.  /  sec.  Omit  if 

gust  is  not  desired. 


Omit  the  gust  card  if  the  gust  is  not  desired.  (Items  53  through  56) 


57  LIlVx  —  Summation  limits  for  chordwise  integration-  Usually 

LIM  1  is  the  first  box  in  a  strip  and  LIM2  is  the  last 
box  where  all  boxes  are  numbered  consecutively. 


58  NOF  —  Panel  number- 


59, 60, 61 
62,  63, 64 
65 


IS1  IS2  IS3 
IS4  IS5  156 
IS7 


AIC  degree -of-freedom  selector.  If  alternative 
#1  is  used  (i.  e.  if  NPC  =  0)  then  IS1  through  IS6 
(IS7  not  used)  represent:  plunging,  pitching, 
control  surface  rotation,  tab  rotation,  control 
surface  plunging,  tab  plunging.respectively.  If 
Alternative  #2  is  used  then  IS1  through  IS7 
represent:  three  cambering  modes  on  the  main 
surface,  control  surface  rotation,  tab  rotation, 
control  surface  plunging, and  tab  plunging, 
respectively. 


66  XHI  —  Reference  axes  measured  in  fractions  of  local  chord 

on  the  inboard  edge  of  a  panel.  (Further  description 
following  XHO-  )  Ignore  if  AIC's  are  not  needed 


67  XHO  ~  Reference  axes  measured  in  fractions  of  local  chord 

on  the  outboard  edge  of  a  panel.  For  Alternative  #1 
5  quantities  are  needed:  the  elastic  axis,  the  leading 
edge  of  the  control  surface,  the  rotation  point  of  the 
control  surface,  the  tab  leading  edge  and  rotation  point, 
respectively.  The  data  pair  XHI  and  XHO  are  repeated 
for  each  panel. 


1* 


Item  No. 


Description 


Omit  data  Items  #58  through  #67  if  AIC's  are  not  required. 


68  NM,  NMD  —  Total  number  of  modes,  i.  e.  ,  max  NQ. 

69  NT,  NTA  —  Total  number  of  ARQ  values. 

70  N5  —  If  "l",save  the  quasi-inverse  on  tape  so  that  future 

modes  may  be  considered.  Otherwise  set  to  "0". 

71  N6  —  Type  of  mode  input  flag 

0  -  polynomial 

1  -  input  h,  dh/dx  from  tape 

2  -  input  h,  dh/dx  from  cards 

72  N7  —  if  "  1 ",  calculate  pressure  forces  and  moments. 

Use  "0"  otherwise. 


7  3  NTPI  -Quasi  Inverse  tape  identification.  Written  on  tape 

when  Quasi  Inverse  is  formed  and  used  to  identify 
the  Quasi  inverse  when  it  is  retrieved  for  additional 
mode  cases. 


74  NT  PM  -  Identification  of  input  mode  (h,  dh/dx)  tape. 


The  input  to  follow  is  associated  with  the  input  modes.  Equations  will  be 
written  to  define  the  variables.  For  panels  (lifting  surfaces)  the  deflections 
will  be  normal  to  the  surface  while  for  bodies  the  deflections  will  be  in  the 
Z-  (if  any)  and  Y-  (if  any)  directions.  For  the  "NFl"^  panel  or  body  the 
deflections  in  the  "NQ"*'  mode  are  calculated  as  follows: 


(f)  (NR> 
NQ 


COEF 


(NR) 


N  -  0,  1 . 


£ 

M  =  0,1,' 


ARQ 


(NR.  NQ. 


where  r  is  in  a  radial  direction  The  origin  of  the  radius  is  either  at  the 
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I 


origin  of  coordinates,  when  N8  =  0,  or  at  the  inboard  edge  (for  the  panel) 
or  the  axis  (for  a  body)  if  N8  is  set  to  1 .  * 

r  =  yj  (y  -  (NSJY^1^)2  +  (z  -  (N8)Zj (NR))2  for  panels 
T  =  ij  (y  -  (N8)Yc(NR))2  +  (z  -  (N8)Zc(NR))2  for  bodies 


\ 


SKETCH  2.  1 .  2-3 


Since  each  panel  is  planar,  '%  is  a  s panwise  distance  in  the  plane 
measured  from  its  inboard  edge  or  from  the  origin  v  ?  -  0:  - 

is  a  radial  distance  only  in  the  sense  that  each  panel  may  have  a 
different  dihedral.  1  he  use  uf  this  radial  distance  is  not  meant  to 
in. ply  that  the  panels  are  curved. 
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Item 

No.  Description 

75  MP,  NMTP  —  Total  number  of  ARQ  coefficients  for  panels  only. 

76  MB,  NMTB  —  Total  number  of  ARQ  coefficients  for  bodies. 

77  NR  —  Panel  number  and/or  body  number. 

78  NQ  —  Mode  number. 

79  NA,  NARQ  —  Number  of  terms  (ARQ)  in  mode  NQ,  panel  and/or  body  NR. 

80  N8  —  Flag  that  sets  origin  of  radial  variable  If  "O' ,  the  x-axis 

is  used.  If"  1 ",  inboard  edge  or  axis  of  body  is  used  as 
origin  (see  equation  for  t). 

The  four  items  NR,  NQ,  NA,  N8  are  repeated  for  each  panel  or 
body  until  all  panels  and  bodies  are  mentioned.  Even  if  a  panel 
or  body  is  not  in  motion  (NA  =  0)  the  set  of  four  items  must  be 
given. 

The  direction  of  motion  of  a  body  is  given  by  the  NZ  or  NY  flags. 

Thus  if  body  pitching  is  given  and  NY  is  flagged,  thei  the  body 
will  be  given  a  yaw.  If  a  body  can  move  in  both  the  Z  and  Y 
directions  (a  nacelle)  then  it  must  be  replaced  with  two  bodies; 
one  with  NZ  =  1  and  one  with  NY  =  1 . 

81  N,  LARQ(l)  —  Power  of  x/s  in  the  mode  polynomial. 

82  M,  LARQ(2)  —  Power  of  rls  in  the  mode  polynomial. 

ARQ  —  Coefficient  of  (x/s)^  (t/s)^  in  the  mode  polynomial. 


Omit  modal  data  cards  if  modes  are  not  desired  (Items  68  through  83). 
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If  modes  are  input  using  the  options  N6=l,  or  N6-2,  omit  Items  73 
through  83,  and  add  the  following  data  items: 


Item 

No, 


84 

NB£X 

Total  number  of  boxes  on  all  panels. 

85 

B 

-  Generalized  force  integration  matrix  (box  area) 
(displacement  1 /4-chord  of  box). 

X 

86 

H_ 

-  Displacement  of  3/4-chord  point  normalized  by 
semispan,  h/s. 

87 

DH1 

-  Derivative  of  normalized  displacement  by  (x/s)  at 
3/4-chord  poin*:  of  bc-.c  .  d(h/s)/d(x/s). 

88 

NBEL 

-  Total  number  of  body  elements. 

89 

B_ 

-  Generalized  force  integration  matrix  (body  element 
length)  X  (diamet  r)  X  (displacement  1/2  -chord  of 
element). 

90 

H 

-  Displacement  of  1  / 2  -chord  point  normalized  by 
semispan,  h/s. 

91 

DH1 

-  Derivative  of  normalized  displacement  by  (x/s) 
I/2-chord  point  of  element.  d(h/s)/d(x/s). 

at 

92 

DH2 

-  d2(h/s)/d(x/s)2. 

NOTE:  If  N6=l,  Items '#84  through  §92  are  input  from  cards  according  to 

the  formats  shown  on  page  8. 

If  N6-2,  Items  §84  through  §92  are  read  from  tape  unit  MTAPE, 
as  one  record  per  variable,  NB£)X,  (B  ,  i  -  1,  NB(DX),  (H  , 
i  =  1,  NBOX)  and  (DH1 i  =  1,  NBOX)l  If  bodies  are  also 
present,  five  more  variables  are  read  from  tape  unit  MTAPE; 
again,  as  one  record  per  variable:  NBEL,  (B  ,  i  =  1 ,  NBEL), 
(EL,  l  -  1,  NBEL),  (DHL,  i  =  1,  NBEL)  and  (L?H2.,  i  =  1,  NBEL). 

For  each  additional  mode,  repeat  Items  585  through  #92. 
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2.  2  Example  Case 


The  following  case  is  to  be  viewed  only  as  an  example  of  the  proper 
input  procedure  and  not  as  an  example  of  an  ideal  or  optimum  idealization. 
The  configuration  is  shown  in  Sketch  2.  2  - 1 . 


c  2.  0 

k  1.0 

r 

M  0.  8 

A  8.  0 

R  6.  0 


Sketch  2.  2-1 

The  fuselage  is  divided  up  into  10  equal  elements  for  the  slender  body 
calculation.  The  divisions  along  with  the  radius  is  given  as  follows: 

F,  =  0,  1.0,  2.0,  3.0,  4  0.  5.0,  h.  0.  7.0.  80.  TO.  10  0 

k 

RAD,  0,  0.7,  0.8,  0  0.  0  95,  1  0,  0.  T  0  8,  0  7.  0  5.  0.0 
k 

The  axis  of  the  fuselage  lies  on  the  x-axis.  7.  V  0  and  it  is  assumed 

that  the  fuselage  moves  only  in  the  /-direction  (N  /’  1).  The  modes  of 
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oscillation  are  symmetric  (6  -  1.0)  pitching  about  the  middle  of  the  fuselage 
and  plunging 

h/s  -  -5.0  +  5  x/8  (where  s  -  5) 

h/  s  =  1.0 

These  modes  apply  to  the  entire  configuration.  The  panels  are  numbered 
from  1  to  5  in  Sketch  2.  2-1. 

The  panels  are  given  as  fallows: 


(1) 

X1 

-  4.  0, 

X2  r 

7.  0, 

X  3  = 

6.  0, 

X4 

=  7 . 

0 

Y1 

=  1.  0, 

^  2 

:  5.0, 

Z1 

0  , 

Z2 

'  --  0 

0 

-  0,1 

0.  333, 

0.  666 , 

1.0, 

T  r 

0,  0 

.  25, 

0.  50,  0.75,  1.0 

(2) 

X1 

=  9  . 

X2 

10  , 

X3  = 

9.  5, 

X4 

10. 

0 

Y. 

=  1.0, 

V2  = 

3.  0, 

Z1  ” 

0  , 

Z2  : 

:  0 

e 

0,  0. 

5,  1.  0 

r 

:  0,  0. 

5.  1.0 

(3) 

xi 

-  2.0, 

X2  - 

10.  0, 

X  3  - 

2.  0, 

X4  = 

10. 

0 

Yi 

=  0  , 

Y  2  = 

1  0, 

Z1  = 

1.0, 

Z2  - 

1  .  0 

e 

-  0,  0. 

'2  5,  0. 

25,  0. 

375,  0. 

50.  0. 

625. 

0.  73. 

0.  873,  0.  ^373. 

0,  1.0 

(4)  X  j  .  X,,  X  ^ ,  X  ,  ,  G  same  as  P) 

Y,  1.0,  Y2  1.0.  /  1.0.  Z,  -10 

r  0,  0.  5.  1  0 

(5 )  X  .  ,  X  , ,  X  ,  X  .  ,  0  ,  r  same  as  i  t ) 

V,  1.0.  V,  0.0.  X  i  -1  .  /  }  -  1 

I  he  input  sheets  ire  shown  on  the  next  several  paces  I  he  output  is  shown 
on  the  subsoquen:  pages 

2 .  2  1  Input  Shoe' s 
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PROGRAM  no 


;.’C  NO"  P'.JNCH  BLANK  COLUMNS  _  NO  l  NDERPUNCH 


PANEL 

MODES 


•iCNS  FOR  XEYPUNCH  00  NOT  PUNCH  ELANK  COLUMNS _ NO  UNDERPUNCHES  IN  SIGN  FIELDS. 


2.  2.  2  Out  put 

#»*C  HIT  dll  tft*tt*tOOVtitt 

••  Matt  0#  atoucm  mimikiii  •• 


1.00*004 


aia.  CMQRO  •  t.MOM  Mf,  tl*l*SP4N  • 

•UCN  HO*  •  o.toooo  at  14  ■  0.40000 


$.00000  iff.  4414  •  ||,t 


•  •• 

4ANIt  NO. 

1  IN4UT  VAIUCS 

*4* 

A 1  • 

4,000000 

*2  • 

1.000000 

41  •  1.000000 

21  • 

0.0 

II  • 

4,000000 

44  • 

4,000000 

42  •  *.000000 

1/  • 

0.0 

NC  • 

4  NS  4  S 

NOttf 

•  l  Nil.  U4 

4 ANILS  •  $ 

4  f»mo»in  oi vis kims  for  muh  i 


0.11»)2**Rf  Oil  0.444444*11  00 

$  *PAN*|*f  DIVISIONS  FOR  444ft 


O.IOOOOOOOf  01 


0.0 

0.2*0000004 

00 

O.IOOOOOOOf  00 

0, 7*0000001 

00 

o.ioooooooe 

01 

20 

•  ••  XI  tUHNIS 

404  4ANEL  !••• 

0.  *00000004 

01 

0.4«***«*5f 

01 

0.4(200001  Of  01 

O.IOOOOOOOf 

01 

0.4*0000001 

01 

0. *11*17**4 

01 

01 

0.100000004 

01 

'*.*00000004  01 

0.54444***4 

01 

0.41*1 I140E 

01 

O.IOOOOOOOf 

01 

0. **0000004 

01 

0 . *******  If 

01 

O.MOOOQOOf  n  1 

O.IOOOOOOOF 

01 

0.400000004 

01 

0 .41)1)7*74 

01 

0.444444 'Of 

01 

O.IOOOOOOOf 

01 

20 

41 A  Elf Nf NTS 

FOR  4ANCL 

O.IOOOOOOOf 

01 

O.IOOOOOOOf 

01 

O.IOOOOOOOf  01 

O.IOOOOOOOf 

01 

0. 70000000C 

01 

0.200000004 

01 

O.IOOOOOOOf 

oi 

O.IOOOOOOOf 

01 

O.IOOOOOOOf  01 

O.IOOOOOOOf 

01 

0.100000904 

01 

O.IOOOOOOOf 

01 

0. *01000004 

01 

0, *00000004 

01 

0.40000000C  01 

0.40000000E 

01 

0.50000000C 

01 

O.IOOOOOOOf 

01 

0.  *00000004 

01 

0. *00000006 

01 

20  Iff  EltNfNIS  *0«  4ANEI  NO.  | 


0.0 

0.4 

( 

0.0 

0.0 

0.0 

0.0 

o.c 

I 

0.0 

o.n 

0.0 

0.0 

o.c 

* 

0.0 

0.0 

0.0 

0.0 

o.c 

1 

4  •••  Xl-J  It INf NTS 

FOR  4ANFL  !••• 

0.42*000001  01 

0.41500000E  01 

0.52*000001 

01 

O.IMOOOOOf 

01 

C. NICOLES  4(1*  4ANEI  NO.  1  ••••• 

0.27500000C  01 

0.22****>*0f  01 

0. 1 7500000C 

01 

0.12*000004 

01 

•  •• 

4ANFL  NO. 

2  INPUT  VALUES 

««• 

XI  •  *.000000 

X2  • 

10.000000 

VI  ■  1.000000 

21  • 

0.0 

11  •  *.*00000 

X*  « 

10.000000 

V2  •  I. 000000 

12  » 

0.0 

NC  •  1  NS  •  1 

NOflf 

■  l  NO.  04 

PANELS  ■  5 

0.0 


0.0 

0.0 


1  CHJROmISE  01  VISIONS  404  FANfL  2 

0.0  0.5OOOOOOOC  00  0. 103099006  01 

!  S4ANMISC  DIVISIONS  404  FANIL  2 


0.0 

O.IOOOOOOOf 

00 

O.IOOOOOOOf  01 

*  •••  XI  fit  Nf NT  S 

FOP  PANEL  !••• 

0. *00900 OOf 

01 

0. *50000004 

01 

O.IOOOOOOOf  02 

0. *25000004  01 

0. *42900001  01 

O.IOOOOOOOf 

02 

0. *50000004 

01 

0. *T500000f 

01 

O.IOOOOOOOf  02 

*  •••  114  ElfNCNTS 

FOR  PANEL  2*** 

O.IOOOOOOOf 

01 

O.IOOOOOOOf 

01 

O.IOOOOOOOf  01 

0.200 00 OOut  01 

0.200000004  01 

0.200000001 

01 

O.IOOOOOOOf 

01 

0. 10000000C 

01 

O.IOOOOOOOf  01 

*  Iff  f LI Nf NTS 

FUR  PANEL  NO.  2 

0.0 


0.0 


0.0 

0.0 


0.0  0.0 
0.0 

2  •••  ll-J  lll*CNT$  404  4 ANCt 


0. *17500001  01  0.414400004  01  - 

C.WiectfS  404  4ANCI  NO.  2 
0.4  4*000004  00  0, *4*000004  00 

•  PANEL  Nil.  *  INPUT  VALUE*  ••• 

«t  *  2.000190  «2  •  lO.OOOOOfl  T|  •  0.0  PI  •  1.000000 


,r>V 
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J.OOOIV)0  It 


II  • 

•K,  «  l<3  »»5 


1tt»U90s»U.J  11  •  1.00JU04  //  •  1.000900 

NUlt  I  •  l  «m |,  uf  »ANt|  *  .  ) 


10  iMiUlMUl  Uimi'JNS  tU«  M<llt  | 


0.0 

o.  Moowwof  oo 


O.l/SOOOOiH  00  0.  »VJOOOOO»  00 

u. •  r^ooooof  jy  o.oiMoooot  oo 

f  SM*rnm  UlvOlUNS  foil  PANU  j 


0.  1  flQOQQOt  00 
0.  lOOOOOO'Jf  01 


u,u 

3,200i>00(Wt  oi 
o.  400110.1001  Oi 
O.IOOOOOOOt  01 

o, •♦loonooot  oi 

0.0 

0.0 

O.lOOOOOOOf  0| 
0. lOOOOOOrtf  01 


o.  iooooooot  0| 
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2.3 


Blank  and  Labeled  Common  Blocks 


2.  3.  1  Blank  Common  Block 


Mnemonics  I  Dimension  Symbol  Source  Description/Equation 


1  NCNSM1, 
NB0X* 


5  N0PAN, 
NP 


1 - 

6 

IQ 

7 

IR 

8 

JSPECS, 

JSP 

9 

FMACH, 

M 

10 

r 

NCARAY 

PART  1  Total  number  of  boxes 
(points)  on  all  lifting 
surfaces  (panels) 

N 

Y  (nc.-l)  (ne.-l), 


where  N  =  N0PAN  = 
number  of  panels 


MAIN  Total  number  of  bodies 
input 


MAIN  Symmetry  flag  -  input 


MAIN  Data  flag  -  input 


MAIN  Total  number  of  panels 
input 


initialized.  Index  of  receiving  point 

in  MAIN, - 

used  in  ,  ,  , 

Part  2  mdex  of  sending  point 


MAIN  Second  symmetry  flag 
input 


MAIN  Mach  Number  -  input 


PART  1  Array  of  the  number  of 
the  chordwise  divisions 
for  all  panels 


A  .lash  through  the  O  indicates  the  letter  O  rather  than  the  number  zero. 
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Item 

No. 


Mnemonics  Dimension 


Symbol 


Source 


Desc  ription/ Equation 


1 1 

NSARAY 

20 

PART  1 

Array  of  the  number  of 
the  spanwise  divisions 
for  all  panels 

12 

N BAR AY 

20 

nba. 

l 

PART  1 

i 

nbai  r  £  <r.c-l)  (ns  -1), 

j--l  J  J 

where  i  =  1,  N0PAN 

13 

ACAP, 

REFARA 

A,  S 

MAIN 

Reference  area  -  input 

14 

B2 , 

S 

b/2,  s 

MAIN 

Reference  semicpan-input 

15 

FL. 

REFCHD 

c.  2b 

r 

MAIN 

Reference  chord  -  input 

16 

BETA 

P 

MAIN 

(1  =  J  1  -  M^  ,  whe  re 

M  -  Mach  Number 

17 

PI 

7T 

MAIN 

rr  =  3.  141  5926, 
built-in  constant 

18 

KR, 

RFREQ 

k 

r 

MAIN 

Reduced  frequency  * 
input  as  RFREQ(20) 

19 

KRDBR 

k  /b 
r  r 

MAIN 

Normalized  k 

r 

20 

GMA 

50 

Vi 

PART  1 

Array  of  dihedral  angles 
for  all  panels  in  radians; 
i  =  1,  N0PAN 

21 

X 

400 

X 

PART  1 

Array  of  x-coordinates 
of  receiving  points 
(3/4  chord  point  of  box)  - 
see  sketch  or.  next  page. 
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Item 

No. 


Mnemonics  Dimension  Symbol 


Source 


Description /Equation 


22 

Y 

— 

400 

y 

PART  1 

Array  of  y~coordinates  of 

receiving  points 

23 

Z  Z 

.  .  .■  - 

400 

. 

z 

PART  1 

Array  of  z-coordinates 
of  receiving  points 

24 

Zl 

400 

Sel 

PART  1 

Array  of 
x -coordinates  ' 

25 

PI 

400 

Vl 

PART  1 

Array  of 
y-coordinates 

1  °f 

inboa  rd 

/  sending 

points 

26 

ZZ1 

400 

V- 1 

PART  1 

Array  of 
z-coordinates 

27 

Z  2 

400 

42 

PART  1 

Array  of  v 

x-coordinates 

28 

P2 

400 

\2 

!^ART  1 

Array  of 
y-coordinates 

^  of 

outboard 
r  sending 
points 

29 

ZZ2 

400 

4-2 

PART  1 

Array  of  , 

7. -coordinates 
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Item 

No. 

Mnemonics 

Dimension 

Symbol 

Source 

Description /Equation 

30 

£V 

400 

iy 

PART  1 

Array  of 

x-coordinates  ^ 

'  of 

31 

PV 

400 

nv 

PART  1 

Array  of 

center- 

y-coordinates  . 

'  line 
sending 

32 

ZV 

400 

4V 

PART  1 

Array  of 

points 

z -coo rdinates  * 

33 

SDELX 

400 

Ax 

s 

PART  1 

Array  of  average 
chordlengths  of  all  boxes 

34 

DELY 

400 

Ay 

PART  1 

Array  of  box  widths. 

1 

See  sketch  above 
items  no.  21-34. 

j 

for 

The  Blank  Common  is  used  by  the  lullowing 
subroutines  of  the  program.  MAIN.  PAK'l  1  ,  GENQ, 
AUGVV  and  PART2. 
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2.3.2  Labeled  Common  Blocks 


Common 

Name 

Source 

User  Subroutines 

/  B0DY / 

PARTI 

GENQ,  AUG W ,  GENF 

/M0DE/ 

GENQ 

AUGW 

/  DLM/ 

KERNEL 

INCR0 

/NTPS/ 

MAIN 

GENB,  GENW,  GENQ, 

PART2  ,  SSLECX,  CXSS, 

AIC ,  AERQ,  FINAL,  GENF 

/  PI  GW  / 

PART  1 

GENB 

/xxz/ 

PART  1 

FINAL,  GENF 

/XY  7./ 

PARI  1 

GENQ,  AER0 

/  YZ  Y/ 

C.ENB 

AER0.  FINAL 

The  following  is  a  detailed  description  of  the  contents  of  each 
labeled  common  block 
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Common 

Name  Mnemonics  Dimension  Symbol  Source  Description/Equation 


PARTI  Body  radius  distribution 


R'o  PARI  1  x-derivative  of  body 
radius  distribution 


PAR  11 


nbe.  -  £  <  nf  -  1  ), 

1  J-l  ■' 

where  nf  is  the  number 

i 

of  body  endpoints  for  the 
i  -th  body:  i  1  ,  NB 


Common 

Name 

n 

Mnrrnomc  s 

Symbol 

Desc  r  ipt  ion 

N«  ( 80 ) 

N8 

Flat:  for  spanwi.se  variables 

; 

NARQ(ftO) 

NA 

Number  of  ARQ  values  for  all 
panels  a  id  modes 

M0DE 

I.ARQfl 20,  2) 

(N.  M) 

Array  of  subscripts  of  all 

ARQ’ s 

ARQ(  1  20 ) 

ARQ 

Array  of  AKQ's 

11(400) 

fh| 

jh]  matrix  \ 

DH1  (400) 

|dh/dx ) 

jdh/dx]  matrix  \  input  for 

0112(400) 

[  d^h  Dlx^  j 

22  I  sPrilal 

[d  h  .■  dx**  ]  ma  rix  /  cases  only 

See  Input  Data,  See  2.  I.  2 

Common 

Name 


Mnemonics 


Symbol 


Source 


Desc  ript  ion /Equation 


DLM 

K10 

k/85 

KERN KP 

Planar 
part  of 

|  the  steady 
contribution 
>  to  the  Kernel 

K20 

K  <9) 
K2 

Nonplanar  ' 

part  of  j 

K  1  RI  1 

RIKj 

R t  (AKj)  Tj 

Unsteady 
,  part  of 

Ke  rnel 

Kirn 

ImKj 

1m  (AKj )  Tj  | 

K2RT2P 

RIK, 

RI  (AK2)  T*  1 

K2II2P 

IrnK^ 

1m  (AK/)  T  * 

K  10  11 

K  ( s )  T 

K1  T1 

K20T2P 

x2,sl 

El 

2 

r 

<•  semi  -width  of 
sending  box 

Common 

Name 


Mnemonics 


Dimension  Symbol  Source 


Description/Equation 


NT  PS 

NTP1 

NTP2 

NTP3 

NTP4 

NT  P'3 

NT  P6 

NTP7 

NT  P8 

NT  P9 

NT  PI  0 

MAIN 

Variable  names  for 
logical  tape  units 
used  in  program 

For  present  assign¬ 
ment  of  units  see 

Sec  2.4  (Logical 

Tape  Units) 

i 

Plow 

f 

CT1 

20 

PARTI 

Array  of  chord- 
lengths  for  inboard 
edge  of  all  panels; 

(>1  :  X2-  X1 

CT2 

20 

‘V 

Array  of  chord- 

lengt  hs  for  out  - 

board  edge  of  all 

panels;  C  X  .  -  X , 

r  O  4  3 

TS 

_ 

50 

avj. 

"  4 

T  1  ')*  1 

where  ~  is  the  >-th 

fractional  spanwtse 

division  of  panel; 

i  1.  NST RIP 

Common 

Name 


Mnemonics 


Dimension 


Symbol  Source 


Description/Equation 


xxz 

X0C 

400 

x/c 

PARTI 

x-location  of  load  points 
as  a  fraction  of  local 
chord 

XIJ 

50 

Xle 

x -coordinate  of  leading 
edge  for  center  of 
strips 

XYZ 

YS 

50 

PARTI 

y -coordinate  of 
centerline  of  strip 

DELYS 

50 

Ay  for  st  rip 

ZS 

50 

z-coordinate  of  center- 
line  of  strip 

DELZS 

50 

Az  for  strip 

KIAMMA 

50 

Dihedral  angle  for 
strip,  in  radians 

CVVIG 

50 

■ 

Average  chord  length 
of  strip 

YXY 

X!  A* 

50 

| 

CENR 

(Alternative  #1  AID 

Elastic  axis  (APT,  *2) 

Location  of  first  j 

cont  rol  point 

1 

X  I  A* 

50 

Location  of  third  control 
point 

X  )A» 

50 

t 

| 

Location  of  rotation 
axis  for  control  surface 

XT  A* 

5  ) 

} 

1 

i 

1 

Location  of  rotation 
axis  fo r  tab 

i 

X2A* 

so 

. . j. 

Location  of  second  j 

control  poir-t  j 

*  All  coordinates  given  in  fractions  of  local  chord 


LOGICAL  TAPE  UNITS 


2.  4 


Tape  Number 

Symbol 

Present 

Assignment 

User  Subroutines 

NT  PI 

1 

MAIN.  GENQ,  PART2,  SSLECX.  S0LVIT, 

QUAS,  FUTS0L 

NTP2 

i  3 

MAIN,  QUAS,  FUTS0L 

NTP3 

3 

MAIN,  GENW,  GENQ,  SSLECX,  CXSS, 

S0LVIT,  QUAS,  FUTS0L,  AIC.  FINAL,  GENF 

NT  P4 

4 

MAIN,  GENQ,  CXSS,  S0LVIT,  QUAS.  FUTS0L 

12 

AUG W, GENF 

NTP8 

8 

MAIN,  GENQ,  QUAS,  FUTS0L 

NT  P9 

9 

MAIN,  GENQ,  S0LVIT,  QUAS,  FUTS0L, 

AIC,  AER0 

NTPIO 

10 

MAIN,  GENB,  GENW,  CENQ,  SSLECX. 

CXSS,  AIC,  FINAL,  GENF 

LTAPE 

1  1 

QUAS,  FUTS0L 

MTAPE 

7 

GENQ 

NAT 

NTP2 

QUAS.  FUTS0L 

RHSTAP 

4 

QUAS,  FUTS0L 

Logical  tape  units  5  and  6  are  used  throughou  the  program  as  the 
standard  READ  and  WRITE  units  respectively. 


Logical  tape  units  MT APE,  NT P9  and  LTAPE  can  be  specified  as 
master  tapes.  MTAPEcan  be  used  as  modal  input  tape  for  special  cases; 
see  2.  5.  6.  4  for  pertinent  information.  NTP9  is  used  for  saving  the  AIC 
supercolumns  when  this  is  desired  —  see  2.  5.  19.  6.  LTAPE  contains  the 
quasi-inverse  matrix  for  future  solutions  of  cases  with  new  modes  only  — 
see  2.  5.  17.6. 


A  \  i 


/n%  % 

u  J6 


iV 
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2.  5  DESCRIPTION  OF  SUBROUTINES 


Program  H7WC  is  operational  on  the  IBM  7094  and  360/65,  and  the 
GE  635  computers.  The  main  part  of  the  program  controls  all  the  twenty- 
one  subroutines  of  the  program:  PARTI,  GENB,  GENW,  GENQ,  AUGW, 
TKER,  PART 2,  INCR&  KERNEL,  IDF1,  IDF2,  SNPDF,  SSLECX.  CXSS, 
SOLVIT,  QUAS,  FUTSOL,  AIC,  AERO,  FINAL  and  GENF.  The  following 
is  a  description  of  each  of  the  subroutines  in  the  order  of  their  use  within 
the  program. 

2.  5.  1  MAIN 


Functional  Description 

The  MAIN  part  of  the  program  reads  the  header  card  for  each  data 
set,  the  Mach  number,  the  reference  area,  the  reference  semi-span 
and  the  reference  chord  as  well  as  the  control  flags  that  define  the  various 
types  of  cases.  MAIN  also  reads  the  array  of  reduced  frequencies,  and 
in  case  of  an  AlC-type  analysis,  additional  control  data  is  read,  as  well 
as  the  gust  data,  if  computation  of  the  harmonic  gust  coefficients  is  also 
desired.  In  case  of  a  modal-type  analysis,  the  control  data  is  read  for 
the  modes  to  be  used.  Depending  on  the  type  of  the  case,  MAIN  selects 
the  appropriate  subroutines  in  an  overall  frequency  -  loop,  except  for 
PARTI  and  GENB  (latter  in  case  of  AIC's  only),  which. are  independent 
of  the  reduced  frequency. 

The  following  is  a  summary  of  all  data  items  used  in  MAIN,  grouped 
according  to  their  source.  Detailed  description  of  each  input  data  item 
is  given  in  Sec.  2.  1.2.  Oth  er  data  items  are  desc  ribed  in  Sec.  2.3 
(Contents  of  Common)  and/or  Sec.  2.4  (Logical  Tape  Units)  as  indicated 
below. 


Data  from  Cards 

FMACH,  ACAP,  FL,  B2,  NDELT,  NP,  NB,  NRF,  NC0RE,  N1.N2, 
N 3,  N4,  RFREQ(20),  NSTRIP,  NB1 ,  NB2,  NPR1  ,  NPR2,  NPR3,  J0BN0, 
J0BN02,  NGUST,  JSPECS,  NPC,  NSV,  NBV ,  NY  AW ,  GZR0,  XZR0, 
VEL,  WG.  LIM{50,  3),  NMD,  NT  A,  N5.  N6,  N7,  NTP1,  NTPM,  NMTP, 
NMTB  • 

See  Sec.  2.  1 . 2  for  dese  ription. 


48 


Data  from  Common 


NCNSM1,  NCARAY (20),  NSARAY{20).  NBARAY(20),  GMA(50), 
X(400),  Y (400),  ZZ(400),  Zl(400),  Pl(400),  ZZ1(400),  Z2(400),  P2(400), 
ZZ2(400),  EV(400),  PV(400),  ZV(400),  SDELX(400),  DELY(400)  - 
from  Blank  Common. 

YS(50),  DELYS(50),  ZS(50),  DELZS(50),  FGAMMA(50),  CWIG(50)  - 
from  Common  /XYZ/. 


XI A(50),  X3A(50),  X5A{50),  X7A(50),  X2A(50)  -  from  Common 
/  YZ  Y/. 

X0C(  400),  CWIG(400)  -  from  Common  /XXZ /. 

CT1  (20),  CT2(20),  TS(50)  -  from  Common  /PI GW/. 

R0  (100),  ROP(IOO),  NBEA(20)  -  from  Common  /B0DY/. 

K10,  K20,  K1RT1,  K1IT1,  K2RT2P,  K2IT2P,  K10IT,  K10T2P, 
E2  —  from  Common  /DLM/. 

See  Sec.  2.  3  for  details. 

Data  from  Tapes 


MNEMONICS 

SYMB0L 

TAPE 

NAME 

DESCRIPTION 

AUGM(530) 

[AWW] 

NTP1 

Matrix  of  complex  downwash 
factors,  read  one  row  at  a  time. 
Dimension  of  matrix  is  NB<^X  x 

NB0X. 

WQ(1  30) 

[Wq] 

NTP3 

Matrix  of  complex  polynomial 
mode  right-hand  sides,  read  one 
row  (of  length  NMD)  at  a  time. 
Dimension  of  matrix  is  NB0X  x  NMD. 

NTP0LD 

LT  APE 

Identification  number  of  case  for 
which  the  quasi-inverse  of  the  matrix 
of  downwash  factors  was  saved  on 
LTAPE  in  a  previous  run.  For 
special  cases  only. 

See  Sec.  2.  4  for  present  assignment  of  logical  tape  units. 


Oat*  to  C ommon 


NB,  NDELT,  NDATA,  N0PAN,  IQ,  IR,  JSPECS,  FMACH,  ACAP, 
B2,  FL,  BETA,  PI,  KR,  KRDBR  —  to  Blank  Common. 

NTPI,  NTP2,  NTP3,  NTP4,  NTP5,  NTP6,  NTP7,  NTP8,  NTP9, 
NTP10  -  to  Common/NTPS/. 

■« 

See  Sec.  2.3  (or  details. 

Data  to  Tapes 


MNEM0NICS 

SYMB0L 

TAPE 

NAME 

DESCRIPTION 

NTPI 

LTAPE 

Identification  number  of  current 
case,  used  to  identify  the  quasi¬ 
inverse  that  is  to  be  saved  on 
LTAPE.  For  special  cases  only, 
when  the  input  data  N5  is  set  to  1. 

AUGM(530) 

[A^WJ 

NTP4 

Augmented  downwash  matrix: 
matrix  of  downwash  factors 
augmented  by  the  matrix  of  poly¬ 
nomial  mode  right  hand  sides 
written  one  row  (of  length  NBOX  + 
NMD)  at  a  time'.  Dimension  of 
matrix  is  NB0X  x  (NB0X  +  NMD). 

2-5.2  SUBROUTINE  PARTI  (JCUM.  CflEFP.  CQEFB,  YIN, 

ZIN) 


Functional  Description 

Subroutine  PARTI  reads  the  geometry  input  of  the  idealized  air¬ 
craft,  the  desired  number  of  spanwise  and  chordwise  divisions  per  panel, 
and  the  fractional  location  of  the  division  lines.  The  panels  are  subdivided 
into  boxes,  and  arrays  of  coordinates  are  generated  for  all  boxes  of  all 
panels.  In  the  case  when  body  interference  is  also  desired,  PARTI  reads 
the  geometry  data  for  all  bodies  and  gcnc-ates  the  necessary  body- 
coordinate  arrays.  The  computed  arrays  are  made  available  for  sub¬ 
sequent  calculations  through  the  Blank  Common  block  placed  in  the 
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appropriate  user  subroutines.  Additional  items  of  computed  geometry 
are  transmitted  Ao  MAIN  via  the  argument  list  of  PARTI,  and  to  several 
other  subroutines  via  the  Labeled  Common  blocks  /P1GW/,  /XX Z/  and 
/XYZ/. 

Data  from  Cards 

XCAP(4),  YCAP(2),  ZCAP(2),  NC,  NS,  C0EFP(2O),  TH(50),  TAU(50), 
ZSC,  YSC,  NF,  NZ,  NY,  C0EFB(2O),  F(50),  RAD(50). 

See  Sec.  2.  1.2  for  details. 

Data  from  Common 

NB,  NDELT,  NDATA,  N0PAN,  IQ,  IR,  JSPECS,  FMACH,  ACAP, 
B2,  FL,  BETA,  PI,  KR,  KRDBR  -  from  Blank  Common  Block. 

See  Sec.  2.  3  for  details. 

Data  from  Tapes 

None 

Data  to  Common 

NCNSM1,  NCARAY{20),  NSARAY(20),  NBARAY (20),  GMA(50), 

X(400),  Y(400),  ZZ(400),  Zl(400),  Pl(400),  ZZ1(400),  Z2(400),  P2(400>, 
ZZ2(400\  EV (400),  PV(400),  ZV(400),  SDELX(400),  DELY(400)  -  to 
Blank  Common  block. 

YS(50),  DELYS(50),  ZS(50),  DELZS(50),  FGAMMA{50),  CWIG(50)  - 
to  Common  /XYZ/. 

X0C  (400),  XIJ  (50)  -  to  Common  /XX  Z/. 

CT1  (20),  CT2(20),  TS(50)  -  to  Common  /P1GW /. 

R0{100),  R0P(100),  NBEA(20),  BGMA(20)  -  to  Common  /B0DY /. 
Data  to  Tapes 

None 
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Pat*  in  Argument  List 


MNEMONICS 

SYMBOL 

DESC  RIPT  ION  /  EQUATION 

JCUM, 

Number  of  strips  on  all  panels 

NSTRIP 

C0EFP(2O) 

Array  of  panel-coefficients  — 
input 

C0EFB(2O) 

Array  of  body-coefficients  —  input 

YIN(20) 

Array  of  y-coordinates  of  inboard 
edge  of  panels 

ZIN(20) 

Array  of  z-coordinates  of  inboard 
edge  of  panels 

2.  5.  3  SUBROUTINE  GENB  (EV.  SDELX.  DELY.  X.  NSTRIP.  BETA, 

IMX.  NPC.  NP.  NSARAY.  LIM) 

Functional  Description 

Subroutine  GENB  reads  the  AlC-type  mode  flags  for  all  panels  of 
the  lifting  surfaces.  It  also  reads  the  array  of  fractional  chord  lengths 
for  the  various  control  points  and  hinge  line  locations  at  the  inboard  edge 
and  at  the  outboard  edge  of  each  panel.  Subroutine  GENB  computes  arrays 
of  coordinates  to  be  used  in  GENB  and  GENW  and  saves  them  on  logical 
tape  unit  NTP10.  Depending  on  the  setting  of  the  flag  NPC  (input  data), 
subroutine  GENB  computes  the  integration  matrix  [B]  according  to  Alter¬ 
native  1  (NPC  =  0)  or  Alternative  2  (NPC  =  1).  (See  Eq.  2.5-12,  Pt  I,  Vol  I.) 
The  [B]  matrix  is  saved  on  logical  tape  unit  NTP10  in  a  partitioned  format, 
and  is  used  in  subroutine  AIC  for  the  computation  of  thcAIC  supercolumns* 
and  harmonic  gust  coefficient  columns,  when  latter  is  also  desired. 

Data  from  Cards 

N0P(2O),  IS(20,  7),  XHI  (20,  ?),  XH0  (20,  7). 

*  A  supercolumn  of  a  matrix  is  a  mat rix -column  whose  elements  are  the 
columns  of  the  original  matrix  taken  in  order;  >.  e.  ,  the  second  column 
is  below  the  first,  the  third  column  is  below  the  second,  et<~  A  square 
matrix  of  order  N  becomes  a  supercolumn  matrix  of  size  x  I. 
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See  Sec.  2.  1. 2  for  detail*. 


Data  from  Common 


NTP1,  NTP2,  NTP3,  NTP4,  NT P5,  NTP6,  NTP7,  NTP8,  NTP9, 
NTP10  —  from  Common  /NTPS/. 

CT1(20),  CT2(20),  TS(50)  —  from  Common  /P1GW/. 

Data  from  Tapes 

None 


Data  to  Common 

XI A(50),  X3A(50),  X5A(50),  X7A(50),  X2A(50)  -  to  Common  /YZY/. 


Data  to  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

NAME 

DESCRIPTION 

NSTRIP 

Total  number  of  strips  on  all 
lifting  surfaces  (panels)-  input  data. 

NSIZB(j.l) 

Number  of  AIC -modes  for 
strip  j 

NSIZB(j,  2) 

NSIZ  B(50,  2) 

Number  of  boxes  in  strip  j 

LIM(50,  3) 

NTP10 

Aerodynamic  box  limits  —  input 

XISLCT(50,  10) 

AIC -type  modal  flags  —  input 

XH(50,  7) 

Array  of  AIC -type  reference 
coordinates 

X(400) 

X 

|  See  Blank  Common 

SDELX(400) 

Ax 

B(30,  10) 

[B] 

Integration  matrix  for  the 
Aerodynamic  Influence 

Coefficients  (AIC's). 
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Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION 

EV(400) 

% 

t 

SDELX(400) 

t> 

X 

a 

>  See  Blank  Common 

DELY  (400) 

Ay 

) 

X(400) 

X 

NSTRIP 

Total  number  of  strips  —  input 

BETA 

p 

See  Blank  Common 

IMX 

Maximum  number  of  rigid  body 
modes  for  case 

NPC 

Camber  flag  —  input 

NP,  N0PAN 

Number  of  panels  —  input 

NSARAY  (20) 

Array  of  the  number  of  spanwise 
divisions  (NS)  for  all  panels 

LIM(50,  3) 

Aerodynamic  box  limits  for  all 
strips  —  input 

2.  5.  4  SUBROUTINE  GENW  (KRDBR.  NCSFTfl.  NPC.  NGUST.  WG. 

qzMIjzMi 

Functional  Description 

Subroutine  GENW  computes  the  AlC-type  substantial  derivative 
matrix,  [W]  described  in  Sec.  2.  5.  1*  .  Depending  on  the  setting  of  the 
flag  NPC  (input  data),  the  [W]  matrix  is  computed  according  to 
Alternative  1  (NPC  =  0)  or  Alternative  2  (NPC  *  1).  The  matrix  is  saved 
for  subsequent  calculations  on  logical  tape  unit  NTP3  in  a  partitioned 
format,  described  in  Sec.  2.  5,  2  *.  In  the  case  when  computation  of  the 
harmonic  gust  coefficient  is  also  desired,  subroutine  GENW  also  compute* 
the  gust  downwash  vector,  and  saves  it  on  logical  tape  unit  NTP3. 

*See  Part  I,  Vol  I 
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Subroutines  GENB  and  GENW  are  used  only  if  an  AlC-type 
analysis  is  desired. 


Data  from  Cards 


None 


Data  from  Common 

YS(50),  DELYS(50),  ZS(50),  DELZS(50),  FGAMMA(50),  CWIG(50)  - 
from  Common  /XYZ  /  . 

Contents  of  Common  /NTPS/. 

Data  from  Tapes 

NSTRIP,  NSIZB(S0,  2),  LIM(50,  3),  XISLCT(50,  10),  XH(50,  7), 
X(400\  SDELX(400)  —  fromNTPlO,  written  in  GENB. 

Data  to  Common 

None 


Data  to  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

NAME 

DESCRIPTION 

GRHS(400) 

!W  ! 

1  gust | 

NTP3 

Gust  right  hand  sides  (downwash 

vector)  W  .  =  cos(r  ->.)• 

kr  R  » 

exp(-ir — (x.-x  )  ),  where 
br  l  o 

T  :  gust  field  declination  from 
*  the  vertical  —  input 

Y.  =  dihedral  angle  of  lifting  sur¬ 
face  at  point  i 

kr  -  reduced  frequency  —  input 

b  reference  sennehord 

r 

x  -  x -coordinate  of  downwash 

i 

point  i 

x  gust  reference  coordinate  — 

o 

input 
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MNEMONICS 

SYMBOL 

TAPE 

NAME 

DESCRIPTION 

W (30,  10) 

tWAIC^ 

NTP3 

Partitioned  AlC.-type  downwash 
matrix.  See  Sec.  2.  5.  2* 

Data  in  Argument  List 


MNEMONICS 

SYMB0L 

DESCRIPTION 

KRDBR 

k  /b 
r  r 

NCSFT0, 

NB0X 

Total  number  of  boxes  on  all 
panels 

NPC 

Camber  flag  \ 

NGUST 

Gust  flag  1 

WG 

W 

g 

Harmonic  gust  amplitude  l  input 

GZR0 

Pg 

Gust  angle  of  declination  I 

XZR0 

X 

0 

Gust  reference  coordinate  / 

•See  Part  I .  Vol  I 
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2.5.  5 


SUBROUTINE  GENQ  (N4,  N5,  N6.  NRF.  JRF.  NBE.  NMD. 

HMI £. JKPflL  MT 

ZIN.  CpEFP.  CQEFB) 

Functional  Description 

Subroutine  GENQ  generates  the  normalwash  boundary  condition  for 
all  lifting  surface  elements.  It  also  generates  the  integration  matrix  used 
to  obtain  the  generalized  forces.  The  normalwash  (r  ranges  over  all 
lifting  surface  box  receiving  points  while  i  ranges  over  all  modes)  is  written 

on  tape  row  by  row  (r  =  1,  i  =  1,2 . NMD,  r  =  2,  i  =  1 , 2 . NMD,  etc.  ) 

to  facilitate  its  use  with  subroutine  S0LVTT.  However,  if  the  number  of 
modes  times  the  number  of  unknowns  (NC0RE)  exceeds  the  size  of  the  work¬ 
ing  array  in  S0LVIT,  and/or  if  the  input  flag  N5  is  set  to  1,  the  W  matrix  is 
written  on  tape  in  column  order,  preceded  by  the  number  of  modes  (NMD), 
ready  for  use  by  subroutines  QUAS  and  FUTS0L.  The  integration  matrix  Br- 
is  written  on  tape  in  column  order. 

The  formula  for  w  used  is  given  by  equations  2.  3-2*and  2.  3-jffor 
cases  without  bodies,  and  by  equation  2.  2-13*for  body  cases;  the  contri¬ 
bution  of  the  body  effect  to  the  normalwash,  wgf,  is  computed  in  sub¬ 
routine  AUGW,  which  is  callea  from  subroutine  GENQ  for  body  cases  only. 
Wgr  is  computed  one  column  at  a  time,  and  is  transmitted  to  GENQ  via 
the  argument  list  of  subroutine  AUGW.  The  total  normalwash  matrix  for 
body  cases,  ( -  w^)  ,  is  then  written  on  tape  in  either  row,  or  column 
order  depending  on  the  size  (NC0RE)  of  the  problem  and/or  on  the  setting 
of  the  input  flag  N5. 

Data  from  Cards 

NR.  NQ,  NARQ(80),  N8(80),  LARQfl  20,  2)  and  ARQ(1 20)  for 
standard  polynomial  cases  (N6  -  0). 

When  modes  are  input  from  cards  (N6  =  2),  the  data  input  from 
cards  are: 

NBOX,  BQ(400),  H(400)  and  :  H1(400). 


*See  Part  I,  Vol  I 
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Data  from  Common 


Contents  of  Blank  Common  and  the  labeled  common  blocks 
/NTPS/  and  /BODY/. 

Data  from  Tapes 

N6  -  0  cases: 

None,  if  NRF  =1.  If  NRF  >  1  and  JRF  >  1  (see  Sec.  2.  5.  5.  ),  the 
following  data  items  are  read  from  tape  NTP10; 

BQ(400),  N8(80),  NARQ(80),  LARQ(120,2),  andARQ(120). 

N6  =  1  cases:  modes  are  input  from  tape  as  follows. 

NB0X.  BQ(400),  H(400)  and  DH1{400)  from  tape  MTAPE;  present 
assignment  is  NTP8  =  8. 

Data  to  Common 


N8(80),  F-'RQ(80),  LARQ(120,  2),  and  ARQ(120), 
H(400),  .11(400)  and  DH2(400)  to  common  /MODE/. 


Data  to  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

BQ(400) 

[Bqi 

NTP10 

Integration  matrix  for  modal 
analysis  written  in  column 
order.  Dimension  of  matru 
is  NT0T  xNMD,  where  NT0T 
is  the  total  number  of  elements 
(-'  ...  NT0T  =  NB0X,  .f  there 
are  no  body  elements,  and 

NT0T  *  NB0X  +  NBE  otherwise- 
where  NBE  -  total  number  of 
body  element s). 

NMD 

NM 

NTP4 

Number  of  modet;  written  only 
•  fM?  1  and -or  NC0RK  5  4000 
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MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

WQ(400) 

fWq1 

NTPW 

Matrix  of  polynomial  mode  right 
hand  sides  (normalwash). 
Dimension  of  matrix  is  NB0X  x 
NMD.  The  matrix  [W  1  is 
written  on  tape  in  row^rder, 
unless  N5  =  1,  and/or 

NC0RE  2  4000;  in  latter  caBe, 

[W  ]  is  written  on  tape  in  column 
orcSer. 

Note  that  NTPW  =  NTP3  in  case 
of  row  order,  and  NTPW  =  NTP4 
in  case  of  column  order. 

BQ(400) 

[Bq’ 

The  following  items  are  written 
on  tape  only  if  more  than  one 
reduced  frequency  is  prescribed 
for  the  case  (NRF>1). 

N8(80) 

NTPB 

f  See  Input  Data,  2.1.2 

NARQ(80)' 

LARQ(120,  2) 

ARQ(120) 

NA 

(N,  M) 

ARQ 

|  Note  that  NTPB  =  NTP2  for  cases 

I  without  bodies,  and  NTPB  =  10 
'  for  body  cases. 

Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION 

N4 

Print  flag  -  input 

N5 

Flag  to  save  quasi-inverse  -  input 

N6 

Modal  input  data  flag  -  input 

NRF 

Number  of  reduced  frequencies  for 
case  -  input 

JRF 

Sequence  number  of  reduced  frequency 
currently  used  in  the  frequency  do-loop 
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SYMBOL 

DESCRIPTION 

NBE 

Total  number  of  body  elements 

NMD 

NM 

Number  of  modes  -  input 

NMTP 

MP 

Total  number  of  ARQ  -  values  for 
all  panels  and  modes  -  input 

NMTB 

MB 

Total  number  of  ARQ  -  values  for 
all  bodies  and  modes  -  input 

NC0RE 

Total  number  of  unknowns  for  case  - 
input 

KD02 

KD/2 

4000  -  built  in  constant  (MAIN) 

YIN(20) 

Y1 

Array  of  y  -cool  _inates  of  inboard 
edge  ox  panels 

ZIN(20) 

Z1 

Array  of  3-coordin«_i.es  of  inboard 
edge  of  panels. 

C0EFP(2O) 

l  coef 

Array  of  modal  coefficients 

for  panels  )  gee 

C0EFB(2O) 

) 

Array  of  modal  coefficients  *  3-3)1 

for  bodies 

2.  5.  6  SUBROUTINE  AUGW  (NMD.  NB0X.  NBE.  NMTB.  N4.  N6. 

NRF.  JRF.  J.  COEFB.  ASUM) 

Functional  Description 

This  subroutine  gives  the  incremental  normalwash,  wgr»  at  all 
lifting  surface  boxes,  r,  due  to  the  axial  singularities  on  all  bodies,  B. 
Slender  body  theory  is  used  to  obtain  the  strength  of  the  axial  singularities. 
The  formulae  used  in  this  subroutine  are  given  by  Equations(2.  2-8)through 
(2.2-ldP  The  incremental  normalwash  values,  wgr«  are  transmitted  to 
subroutine  GENQ  via  the  argument  list  under  the  mnemonic  name  ASUM 
one  columr.  (for  one  moue)  at  a  time.  This  subroutine  also  computes  the 
inte- ration  matrix  of  all  body  elements  for  the  calculation  of  the  generalized 
fo  s  -ji'  body  cases. 

Subroutine  AUGW  is  bypassed  for  data  without  body  elements, 


*S«e  Part  I,  Vol  I 
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Data  from  Cards 


None  for  the  special  case  when  N6  =  1  (data  input  from  tape  MTAPE). 

For  standard  polynomial  cases  (N6  =  0)  data  input  from  cards  is 
same  as  for  subroutine  GENQ,  except  that  data  is  now  read  for  body  elements 
instead  of  panel  boxes. 

Data  from  Common 

Contents  of  Blank  Common,  and  the  labeled  common  blocks  /B0DY/ 
and  /M0DE/. 


Data  from  Tapes 


None  for  standard  polynomial  mode  cases.  For  the  special  case 
when  the  modal  data  flag  N6  =  1,  the  following  data  items  are  input  from 
tape: 


i 

MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

NBE 

Number  of  body  elements 

BB(80) 

[B  1 

Integration  matrix  for  the 

NTP8 

body  elements 

H(400) 

[h] 

[  h  ]  matrix 

DH1  (400) 

[dh/dx] 

[dh/dx  1  matrix 

DH2(400) 

[d2h/dx2] 

2  2 

[d~h  £x  ]  matrix 

The  above  three  matrices 
are  read  in  column  order. 
Dimension  of  matrices  is 
NBE  x  NMD. 

Data  to  Common 


MNEM0NICS 

SYMB0L 

TAPE 

DESCRIPTION 

BB(80) 

B 

NTP8 

One  column  of  the  integration 

4 

matrix  for  body  elements 
written  on  tape  for  one  mode; 
length  of  column  is  NBE. 

The  following  items  are 
written  on  tape  only  if  more 
than  one  reduced  frequency 
is  prescribed  for  case 
(NRF  >1) 

N8(80) 

N8 

NTPB 

1  See  Input  Data  2.1. 2 

NARQ(80) 

NA 

LARQ(120,2) 

(N,  M) 

( 

ARQ(120,  2) 

ARQ 

/ 

Data  in  Argument  List 


MNEM0NICS 

SYMB0L 

DESC  R1PTI0N 

NMD 

NM 

Number  of  modes  input 

NB0X 

Total  number  of  boxes  on  all 
lifting  surfaces 

NBE 

Total  number  of  body  elements 

NMTB 

MB 

Total  number  of  ARQ  values 
for  all  bodies  and  modes  •  input 
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MNEMONICS 

SYMBOL 

DESCRIPTION 

N4 

Print  flag  for  writing  the 
[APZ],  [APy]  matrices 

N6 

Modal  input  flag 

>  input 

NRF 

Number  of  reduced 
frequencies  for  case  ' 

JRF 

Sequence  number  of  current 
frequency  in  the  frequency 
do-loop 

J 

Sequence  number  of  current 
mode  in  modal  do-loop 

C0EFB(2O) 

coef 

Array  of  scale  factors  fi  r 
body  deflection  modes. 

ASUM(400) 

Incremental  [W]  due  to  body 
interference 

2.5.7  SUBROUTINE  TKER  (XO  .  YQ  .  ZO  .  KR.  BR.  GAMS. 

GAMSIG.  KKR.  KKI,  FMACH) 

Functional  Description 

Subroutine  TKER  evaluates  the  total  unsteady  kernel  for  one 
receiving-point/sending-point  combination  for  body  cases  only.  (See  Eq.  A.  5.)* 
Subroutine  TKER  is  essentially  the  same  as  subroutine  KERNEL,  except 
that  the  total,  rather  than  the  incremental  oscillatory  kernel  is  computed 
here,  and  that  the  planar  and  nonplanar  parts  of  the  kernel  are  not  separated. 
Subroutine  TKER  is  called  from  subroutine  AUGW;  it  is  bypassed  in  the  case 
when  the  input  data  contains  no  bodies. 


Subroutine  TKER  uses  no  card-,  tape-,  or  common-data  input  or 
output;  data  transmission  to  and  from  TKER  is  done  through  its  argument 
list. 

*See  Part  I,  Vol  I 
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Data  in  Argument  LiBt 


MNEMONICS 

SYMBOL 

DESCRIPTION 

xo 

X 

o 

x.  -  | 

1  J  1 

Y0 

yo 

^i  ’Ij  /  See  Blank  Common 

ZO 

z 

o 

N 

»-»• 

1 

KR 

k 

r 

Reduced  frequency 

BR 

b 

r 

Reference  semi-chord 

GAMS 

Ys 

Dihedral  angle  of  sending  point  j 

GAMSIG 

Dihedral  angle  of  receiving  point  i 

FMACH, 

M 

M 

Mach  Number 

E2 

2 

e  ,  where  e  is  semi-width  of 
sending  box 

KKR 

Real  part  of  total  (steady  + 
unsteady)  kernel 

KK1 

Imaginary  part  of  total  kernel 

2.5.8  SUBROUTINE  PART2  (NYAW.  NBV) 


Functional  Desc  ription 

Subroutine  PART2  is  used  for  the  calculation  of  the  matrix  of 
downwash  factors,  [D],  (See  Eq.  2.  1-14.*)  It  prepares  the  arguments  to 
subroutines  SNPDF  and  INCR0,  which  do  the  actual  computation  of  the 
downwash  factor  elements,  and  calls  these  in  a  do-loop  for  all  boxes  of 
all  lifting  surfaces.  The  steady  part  of  the  downwash  factor  elements  is 
computed  in  subroutine  SNPDF,  and  the  incremental  oscillatory  part  in 
subroutine  INCR0;  the  call  to  INCR0  is  bypassed  for  steady  cases  (kr  =  0.  0). 
Subroutine  PART  2  also  accounts  for  the  effect  of  the  (two)  symmetry  plane(s). 

Data  from  Cards 


None 


*See  Part  I,  Vol  I 
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Data  from  Common 


Contents  of  Blank  Common,  Common  /NTPS/,  and  Common 

/DLM/. 


Data  from  Tapes 


None 


Data  to  Common 


None 


Data  to  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

AWW  (400) 

AWWI(400) 

!  [A-], 

)  ED) 

NTP1 

Matrix  of  downwash  factors 
written  on  tape  in  row  order. 
Dimension  of  matrix  is 

NB0X  x  NB0X. 

The  real  parts  of  the  matrix 
elements  are  called  AWW  in 
the  program,  and  the  imaginary 
parts  are  called  AWWI. 

Data  in  Argument  List 


NYAW 

NBV 


yaw  flag 

number  of  boxes  on  vertical  panels 
(  fin  )  of  T -tail  case 


input  in  MAIN 


2.5.9  SUBROUTINE  INCRO  (AX.  AY,  AZ ,  AX1 ,  AY1 ,  AZ 1 .  AX2,  AY2, 

AZ2,  GAMS,  GAMSIG.  LHS.  IR.  lO.  NBXS. 
NCPNB.  NDBLE,  NBV.  DELR.  PELT.  FL, 
BETA.  SDELX,  DELY.  KR) 


Functional  Description 


Subroutine  INCR0  prepares  the  argument  for  the  subroutines 
KERNEL,  IDF1  and  IDF2.  It  calls  subroutine  KERNEL  to  compute  the 
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incremental  oscillatory  part  of  the  kernel,  K,  for  each 
receiving -sending  box  combination  at  the  three  points  of  the  bound  vortex 
segment:  at  the  center  (Kc),  at  the  inboard  point  (IC)  and  at  the  outboard 
point  (K0).  Since  even  a  relatively  small  case  requires  many  kernel 
computations  (e.  g.  ,  a  wing  with  100  boxes,  in  symmetry,  requires 
100x300x2  =  60,  000  kernel  values),  extra  programming  effort  was  made 
to  save  on  computing  time.  Neighboring  strips  within  a  wing  panel  have 
a  common  kernel  on  the  common  boundary;  this  fact  is  utilized  in  INCR0 
so  that  the  computing  time  of  the  kernels  is  reduced  by  about  30  percent. 

After  the  triplet  of  kernels,  (Kc>  K.,  K0),  is  obtained  for  one 
receiving-sending  box  combination,  subroutine  INCR0  computes  the 
coefficients  of  the  parabolas  for  the  numerical  integrations  done  in  two 
parts.  The  planar  part  of  the  integration  is  done  in  subroutine  IDF1, 
and  the  nonplanar  part  in  IDF2. 


Subroutine  INCR0  uses  no  card-,  or  tape  input  or  output; 
transmission  of  data  is  done  through  its  argument  list  and  the  labeled 
common  /DLM/. 


Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION/ EQUATION 

AX 

X 

o 

x  -  i 

v  V 

AY 

y  -  n  \ 

7  V  1 

AZ 

z 

o 

Z  ;v  I 

AX1 

x  •  &el  f 

AY1 

V  See  Sec  2.3.1 

y  '  nel  / 

{  (Blank  Common) 

AZ  1 

*  •  ^el  y 

AX2 

x  -  ieZ  | 

AY2 

y  *  ,5c2  1 

AZ2 

7  '  ^c2 

/ 
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MNEMONICS 


SYMBOL 


•DESCRIPTION/ EQUATION 


GAMS 

GAMSIG 

LHS 

IR 

10 

NBXS 

NCPNB 

NDBLE, 
JSPECSj  JSP 

NBV 

DELR 

DELI 

FL. 

REFCHD 

BETA 

SDELX 

DELY 


Dihedral  angle  of  receiving  point  i 

Dihedral  angle  of  sending  point  j 

Flag;  0  when  effect  of  right  wing  is 
computed,  1  for  left-wing 

Index  of  sending  point 

An  index  running  from  1  through 
NCM1,  where  NCM1  is  the  number 
of  chordwise  boxes  of  the  strip  in 
which  the  sending  point  lies. 

The  number  of  boxes  on  the  wing 
panel  in  which  the  sending  point 
lies  +  the  total  number  of  boxes  of 
the  preceding  panels. 

The  number  of  boxes  in  the  first 
strip  of  the  panel  in  which  the  sending 
point  lies  +  the  total  number  of  boxes 
of  the  preceding  panels. 

Second  symmetry  flag  -  input 


Total  number  of  boxes  on  vertical 
panels  of  T-tail  cases  -  input 

Real  part  of  the  incremental 
oscillatory  downwash  factor 

Imaginary  part  of  the  incremental 
oscillatory  downwash  factor 

Reference  chord  -  input 


Average  box -chord  of  sending  box 
2c,  box  width  of  sending  point 
Reduced  frequency 


KR 


Functional  Description 

Subroutine  KERNEL  evaluates  the  modified  ire  remental  oscillatory 
kernel,  FT,  for  one  receiving-sending  point  combination  at  a  time  The 
planar  and  nonplanar  parts  of  the  kernel  are  evaluated  oeparately  and  are 
returned  to  the  calling  program,  subroutine  INCR0,  via  the  labeled 
common  / DLM/.  (Evaluates  right  hand  side  of  Eq.  (2.  1-13).*) 


Subroutine  KERNEL  uses  no  card  or  tape  input  or  output; 
transmission  of  data  is  done  through  the  argument  list  and  the  labeled 
common  /  DLM/. 


Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION/ EQUATION 

XO 

X 

AX  for  center 

o 

AX!  for  inboard  edge  ' 
AX2  for  outboard  edge 

YO 

yo 

AY  for  center 

AY1  for  inboard  edge 
AY2  for  outboard  edge  j 

of  sending 
*  box;  see 

2.  5.  9 

ZO 

z 

AZ  for  center 

1 

o 

AZ1  for  inboard  edge 
AZ2  for  outboard  edge 

KR 

k 

r 

Reduced  frequency 

BR 

b  .C/2 
r 

Reference  semichord 

GAMS 

Yr 

) 

GAMSIG 

Ys 

|  see  2.  5.  9 

MNEMONICS 

SYMBOL 

DESCRIPTION/ EQUATION 

EPS 

€ 

e  =  0.  00001;  built  in  constant 
from  INCR0 

T1 

T! 

T1  =  cos  <vr  '  V 

T2 

rn 

a2 

1 

T2  =  [*oC08Yr -yosinYr  )(ZoCOSYs'yoSinVl/r 

M, 

M 

Mach  number 

FMACH 

2.5.11  SUBROUTINE  IDF1JEE.  E2.  AT 2.  ETA01 .  ZET01 .  ARE. 

AIM.  BRE.  BIM.  CRE.  CIM.  R1SQX. 
XIIJR,  XIIJI) 


Functional  Description 

Subroutine  IDF1  performs  the  integration  of  the  planar  parts  of 
the  incremental  oscillatory  kernels;  see  Eq.  B.  9*  The  result 
of  the  integration  is  the  complex  number  (XIIJR,  XIIJI),  returned  to  sub 
routine  INCR0  via  the  argument  list  of  IDF1. 


Subroutine  IDF1  uses  no  card-,  tape-  or  common  input  or  output; 
data  transmission  to  and  from  IDF1  is  done  through  its  argument  list. 

Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION 

EE 

e, Ay  / 2 

1  s 

2 

Semiwidth  of  sending  box  j 

E2 

e 

AT  2 

|T2| 

See  2.  5.  10 

ETA01 

y 

y  cosY  +  z  sinY 

7o  r  o  r 

ZET01 

z 

z  c o 8  V  -  y  sin  V 
o  r  7o  r 

*See  Part  I,  y0l  I 


MNEMONICS 

SYMBOL 

DESCRIPTION 

ARE 

RefAj) 

N 

AIM 

Im(Aj ) 

1 

BRE 

Re(Bj) 

Coefficients  of  the  parabola 

f  for  planar  part  of 

BIM 

Im(Bj ) 

| 

kernel-integration 

CRE 

Re(Cj) 

1 

CIM 

lm(Cj) 

R1SQX 

2 

r 

-2  A  -2 
y  +  z 

XII  JR 

Real  part  of 
integral  1 

planar 

XIIJI 

Imaginary  part  1 
of  integral 

contribution 

2-5.12  SUBROUTINE  IDF2  (EE.  E2.  AT2,  ETA01 ,  ZETOI .  A2R, 

A2I.  B2R.  B2I.C2R.C2I.  R1SQX. 
PUJR.  PUJI) 

Functional  Description 

Subroutine  IDF2  performs  the  integration  of  the  nonplanar  parts 
of  the  incremental  oscillatory  kernel;  see  Eqs.  (B.  15)  and  (B.  16).* 

The  result  of  the  integration  is  the  complex  number  (DIIJR,  DIIJJ), 
returned  to  subroutine  INCR0  through  the  argument  list  of  IDF2. 


Subroutine  1DF2  uses  no  card  .  tape  or  common  input  or  output; 
data  transmission  to  and  from  IDF2  is  done  through  its  argument  list. 


•See  Part  I,  Vol  I 
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Data  in  Argument  List 


|  MNEMONICS 

SYMBOL 

DESCRIPTION 

EE 

• 

e 

2 

E2 

e 

See  2. 5. 11 

AT  2 

It2I 

i 

ETA01  • 

y 

1 

ZET01 

Z 

) 

A2R 

Re(A2) 

1 

A2I 

Im(A9) 

Coefficients  of  the  parabola 

for  nonpianar  part  of 

B2R 

Re(B2) 

^  kernel  integration 

B2I 

Im(B,) 

La 

f 

C2R 

Re(C2) 

1 

C2I 

Im(C2) 

; 

R1SQX 

See  2.  5.  1 1 

DIIJR 

Real  part  of 

integral  j  nonpianar 

DIUI 

«  .  (  contribution 

Imaginary  part  J 

of  integral 

_ 

2-5.i3  SUBROUTINE  SNPDF  (SL.  CL.  TL.  SGS.  CGS.  SGR.  CGR. 

XO  .  YQ  .  ZO  .  EE.  DIJ.  BETA.  CV) 


Functional  Description 

Subroutine  SNPDF  computes  the  steady  downwash  factors  for  one 
receiving-sending  box  combination  at  a  time.  The  result,  DIJ,  is 
returned  to  PART2  via  it*  argument  list.  (See  Eq.  C  35b,  Part  1,  Vol  1.) 


Subroutine  SNPDF  uses  no  card.  tape,  or  common  data  input  or 
output;  data  transmission  to  and  from  SNPDF  is  done  through  its  argument 
list 
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Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION 

SL 

sin  A,  whereAiB  the  sweep  angle 
of  the  l/4-chord  line  of  sending 

box  [sweep  angle  of  bound  vortex] 
modified  by  Prandtl-Glauertrule. 

CL 

cos  A 

TL 

tan  A 

SGS 

sin  Y 
s 

CGS 

cos  y 

8 

SGR 

sin  Yr 

CGR 

cos  Yr 

XO 

) 

YO 

/  See  2.  5.  10 

ZO 

) 

EE 

e 

Semi-width  of  sending  box 

D1J 

Steady  contribution  to  downwash 
factor 

BETA 

P 

f  7 

•Jl  -  M 

CV 

c  ,AX 

avg  s 

Ax,  average  chord  length  of 
sending  box 

2.  5.  14 


SUBROUTINE  SSLECX  (NCNSM1 .  MD.  NGUST1 
Functional  Description 


Subroutine  SSLECX  is  used  with  AlC-type  data,  and  only  if 
NC0RE  <  4000.  It  reads  the  matrix  of  downwash  factors  [D]  from  tape 
NTP1  one  row  at  a  time,  and  augments  it  with  the  AlC-type  downwash 
matrix  read  from  tape  NTP3.  For  gust  cases,  when  NGUST  #0, 

subroutine  SSLECX  also  reads  the  gust  downwash  vector  from  tape  NTP3. 
The  augmented  matrix  [DjW]  is  saved  on  tape  NTP4  in  row  order,  ready 
to  be  used  by  subroutine  S0LVIT  to  obtain  the  solutions  to  the  system  of 
simultaneous  linear  equations  represented  by  the  matrix  [  Dj  W|. 

Data  from  Cards 


None 


Data  from  Common 
Contents  of  Common  /NT PS/. 


Data  from  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

NSTRIP 

) 

NSIZ  B(50,  2) 

NTP10 

Isee  2.  5.  3 

LIM(50,  3) 

) 

AUGM(530) 

[A] 

NTP1 

Matrix  of  downwash  factors 
written  on  tape  in  row  order 
in  subroutine  PART2.  Dimen¬ 
sion  of  matrix  is  NB0XxNB0X. 

GUST(400) 

{Wgust} 

NTP3 

Gust  downwash  vector  \ 

(  See 

W(30,  10) 

twUic 

NTP3 

AIC  type  downwash  1  2.  5.4 

(normalwash). 
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Data  to  Common 


None 


Data  to  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

AUGM(530) 

[A  j W] 

NTP4 

Augmented  downwash 
factor  matrix 

Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION 

NCNSM1 

n 

Total  number  of  boxes  on  all 
lifting  surfaces  (panels) 

MD 

m 

Total  number  of  AlC-type 

right-hand  sides  NSTRIP 

V  XISLCT. 
m  =  jmax 

j=  1 

where  NSTRIP  =  total  number  of 

strips  on  all  lifting  surfaces,  and 

XISLCT.  =  the  number  of 

jmax 

elastic  modes  for  j-th  strip  - 
whenever  NGUST  =  0.  When  NGUST±0 
NSTRIP 

m  =  2-i  XISLCT.  +  1. 

.  ,  imax 

J  =  1  J 

NGUST 

Gust  flag  -  input 

Functional  Description 


Subroutine  CXSS  is  used  with  AlC-type  data,  and  only  if 
NC0RE  2  4000,  and/or  N5  =  1,  It  reads  the  AlC-type  downwash  matrix, 
[W]^.  from  tape  NTP3  and  forms  columns  of  right-hand  sides.  When 
NGUST  #  0,  subroutine  CXSS  also  reads  the  gust  downwash  vector  from 
tape  NTP3;  this  becomes  the  last  column  of  right  hand  sides.  These  are 
saved  on  tape  NTP4  preceded  by  the  number  of  right  hand  sides  for  the 
case,  sdy  to  be  used  by  subroutine  QUAS  and  FUTS0:L  to  obtain  the 
solutions. 


Data  from  Cards 

None 

Data  from  Common 
Contents  of  Common  /NTPS/. 
Data  from  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

NSTR1P 

) 

NSIZB(50,  2) 

NTP10 

>  See  2.  5. 3 

LIM(50,  3) 

) 

GUST(400) 

Kl 

NTP3 

Gust  downwash  vector 

W  (30,  10) 

Maic 

NTP3 

J  see 

Partitioned  matrix  (2.  5.  4 

of  AlC-type  down-  ' 

wash  (normalwash) 

Data  to  Common 


None 
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' 

MNEMONICS 

SYMBOL 

TAPE 

_ 

- -  - -  - 

DESCRIPTION 

NMD 

m 

>NTP4 

Number  of  AlC-type  right 
hand  sides  (elastic  modes), 
including  the  gust  right 
hand  side  whenever 

NGUST  *  0 

WX(400) 

J 

Matrix  of  AlC-type  right 
hand  sides  written  in 
column  order.  Dimensions 
of  matrix  is  NB0X  x  NMD 

Data  in  Argument  List 


Same  as  in  2,  5,  14 

2.5.16  SUBROUTINE  S0LV1T  (ND.MD.KD.Nl.NM.NGi.NW.NPRn 

Functional  Description 

Subroutine  S0LVIT  is  used  for  cases  with  NC0RE  <  4000, 
and/or  N5  =  0.  It  solves  the  system  of  simultaneous  linear  equations 
represented  by  the  augmented  matrix  [D  j  W],  written  on  tape  NTP4.  The 
solutions  are  saved  on  tape  NTP3  in  column  order;  the  data  on  the  input 
tape,  NTP4,  is  not  preserved. 

Data  from  Cards 


None 


Data  from  Common 


None 


Data  from  T apes 

AUGM(530),  the  augmented  downwash  factor  matrix  from  tape 
NTP4  -  see  2.  5.  1  or  2.  5.  14. 
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Data  to  Common 


None 


Data  to  Tape 8 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

A(4000) 

NTP3 

Solution  matrix  written  in 
column  order;  dimension 
of  matrix  is  NDxMD, 
where  ND  =  NB0X,  and 

MD  =  NMD. 

Note  that  subroutine  S0LVIT  uses  a  single  complex  work  array,  A(4000), 
for  all  tape  input  and  output  variables. 


Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION 

. 

ND, 

NCNSM1, 

NB0X 

n 

Total  number  of  boxes  on  all  lifting 
surfaces  (panels)  =  size  of  the 
square  matrix  [Aww]  =  [D]  =  matrix 
of  downwash  factors. 

MD 

m 

Total  number  of  right  hand  sides  in 
the  system  of  simultaneous  linear 
equations  solved  in  subroutine  S0LVIT. 

KD 

Working  array  size  (real-variable 
dimension);  presently  set  to  8000 

NI 

NTP4;  tape  number  assigned  to 
tape  containing  all  rows  of  the 
augmented  matrix  [Aww  1  W]  - 
see  2,  5.  1  or  2.  5.  1 4 

NM 

NTP1;  tape  number,  used  as 
scratch  unit 

N0 

NTP9;  tape  number,  used  as 
scratch  unit 

NW 

NTP3;  tape  number  assigned  to  tape 
containing  all  solutions  of  the  system 
of  simultaneous  linear  equations 

NPR1 

Print  flag  for  solutions  -  input 
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For  present  assignment  of  tape  units  see  2.4. 


2.5.17  SUBROUTINE  QUAS  (ND.  MD.  KD.  NI.  MM.  N0.  NAT.  NW. 

LTAPE.fiHSTAEvMPRl) 

Functional  Description 

Subroutine  QUAS  is  used  in  conjunction  with  subroutine  FUTS0L 
for  cases  with  NC0RE  £  4000,  and/or  N5  -  1.  It  reads  the  matrix  of 
downwash  factors,  [D]  from  tape  NT  PI  in  row  order,  and  as  many 
columns  of  right  hand  sides  (m.say),  from  tape  RHSTAP,  as  the 
dimension  of  the  work  area  in  QUAS  (3700  complex  words)  will  permit. 
Subroutine  QUAS  computes  the  so  called  quasi-inverse  of  the  matrix  [D] 
and  saves  it  on  tape  LTAPE;  it  also  computes  the  solutions  for  the  m 
right  hand  sides  and  saves  these  on  tape  NTP3. 

Data  from  Cards 


None 


Data  from  Common 


None 


Data  from  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

A(3700) 

[D], 

[AWW] 

NTP1 

Matrix  of  downwash  factors; 
see  2.  5.  8 

NMD 

NTP4 

Total  number  of  modes  for  case 

A(3700) 

[W] 

NTP4 

Matrix  of  right  hand  sides; 
see  2.  5.  5  or  2.  5.  1  5 

Data  to  Common 


None 


Data  to  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

A(3700) 

LTAPE 

Quasi-inverse  matrix  of 
[D].  Master  tape  has  to  be 
specified  to  preserve  quasi¬ 
inverse  for  future  analysis; 
otherwise  tape  is  considered 
a  scratch  unit. 

A(3700) 

NTP3 

Solution  matrix;  dimension 
NB0Xxn,  where  n  is  the 
number  of  right  hand  sides 
for  which  solutions  were 
obtained. 

Note  that  Subroutine  QUAS  uses  a  single  complex  work  array,  A(3700) 
for  all  input  and  output  variables. 


Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION 

ND, 

NCNSM1 , 

NB0X 

n 

Total  number  of  boxes  on  all  lifting 
surfaces 

MD 

m 

On  output  MD  =  n  =  the  number  of 
right  hand  sides  for  which  solutions 
were  obtained 

KD 

. 

Working  array  size  (complex 
variable  dimension);  presently 
set  to  3700. 

NI 

NTP1;  tape  number  assigned  to 
tape  containing  all  rows  of  matrix 
[Aww]  (matrix  of  downwash 
factors) 

MM 

NTP8 

\ 

N0 

NTp9  (  tape  numbers  for 

1 '  |  sc ratch  units 

NAT 

NTP2 
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MNEMONICS 

SYMBOL 

DESCRIPTION 

NW 

NTP3,  tape  number  for  the 
solution  matrix 

LTAPE 

Tape  number  assigned  to  the  tape 
containing  the  quasi-inverse  of 
[D] 

RHSTAP 

NTP4;  tape  number  assigned  to 
tape  containing  the  matrix  ot 
right  hand  sides 

For  present  assignment  of  tape  units  see  Sec.  2.4 


2.5.18  SUBROUTINE  FUTS0L  (ND.  MD.  KD.  NL  MM.  NO.  NAT. 

NW.  LTAPE.  RHSTAP) 


Functional  Description 

Subroutine  FUTS0L  is  used  for  cases  with  NC0RE  a  4000,  and 
whenever  new  modes  only  are  required  for  a  previously  computed  case 
(N1  =  1).  When  N1  =  0,  subroutine  FUTS0L  is  called  following  sub¬ 
routine  QUAS,  and  uses  the  quasi-inverse  of  the  matrix  [D]  from  tape 
LTAPE.  Additional  right  hand  sides  are  read  from  tape  RHSTAP,  and 
the  additional  solutions  are  saved  on  tape  NTP3  following  the  previously 
obtained  solutions. 

Subroutine  FUTS0L  is  bypassed  when  N5  =  1  and  NC0RE  <  3700. 

Data  from  Cards 


None 


Data  from  Common 


None 
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Data  from  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

A(3700) 

[W] 

RHSTAP 

Additional  columns  of 
right  hand  sides 

A(3700) 

LTAPE 

The  quaBi-inverse  matrix 
generated  in  QUAS 

Data  to  Common 


None 


Data  to  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

A(3700) 

NW 

The  additional  solutions  for 
the  additional  right  hand  sides 

For  present  assignment  of  tape  units  see  Sec.  2.  4. 

Data  in  Argument  List 

Same  as  for  QUAS  -  see  2.  5.  17,  except  that  NI  is  now  merely 
a  scratch  unit  and  does  not  contain  the  input  matrix  [D]. 

2.5.19  SUB&flmiNE  1C0P. 

VEL.  NGUST.NPR2) 

Functional  Description 

Subroutine  AIC  is  used  for  AlC-type  cases  only.  It  reads  the 
solutions  to  the  AIC. -type  elastic  modes  from  tape  NTP3  and  the  partioned 
integration  matrix,  fB],  from  tape  NTP10.  Using  these  arrays  subroutine 
•\IC  supercolumns  according  to  Eq.  2.  5-1*.  and  saves  them  on  tape  NTPq. 


*See  Part  I.  Vol  I 
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Data  from  Cards 


None 

Data  from  Comnr 
Contents  of  the  labeled  < 
Data  from  Tapes 


MNEMONICS 

SYMBOL 

NSTRIP 

) 

NSIZE(50,  2) 

LIM(50,  3) 

I 

B(10,  30) 

fBUic 

DUM1 W (400,  4) 

(Dr'(w]AIC 

CPGUST(400) 

Data  to  Common 


None 

Data  to  Tapes 


MNEMONICS 

SYMBOL 

CH(4.  4.  50) 

blocks  /NT PS/,  /XYZ/  and  / YZY/ 


TAPE 

DESC  RIPTION 

NTP10 

See  2.  5.  3 

Integration  matrix 

NTP3 

Solution  matrix  for  the  A1C 
type  right  hand  sides 

NT  P3 

Solution  for  the  gust  right 
hand  6ide  -  read  only  if 
NGUST  *  0. 

TAPE 

DESCRIPTION 

NT  P9 

AIC  matrix  written  one 

supcrcolumn  at  a  time; 

see  Eq  2  b-2.  * 

Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION 

BR 

br>  c/2 

Reference  semichord 

s, 

REFSPN 

S 

Reference  semispan  -  input 

KR 

v2* 

Reduced  frequency 

NBF 

First  bay  number  1  .  ,  .  , 

'  (for  which 

NBL 

Last  bay  number  I  ^*9  s  are 

written  on 

J0BN0 

AlC-case-identifier  taPe 

IC0D 

1 .  built  in  constant 

VEL 

Aircraft  velocity  1 

NGUST 

Gust  flag  \  -  input 

NPR2 

Print  flag  for  AlC's  J 

2.  5.  20  SUBROUTINE  AER0  iBR.  B2.  KR.  ACAP.  NRF.  NSV.  XYAW. 

JMX.NPQ.NDELI.NPR}} 

Functional  Description 


Subroutine  AER0  is  used  for  AlC-type  cases  only.  It  reads  the 
AIC  supercolumns  from  tape  NTP9  and  computes  the  stability  derivatives 
for  the  rigid  body  modes  according  to  Eqs.  (2-4-b)  -  (2.4-121*  This 
subroutine  is  primarily  used  fcr  the  monitoring  of  the  AlC's.  It  can  be 
bypassed  by  setting  the  data  flag  NPR3  to  0. 

Data  f  ror>i  Cards 


None 


Data  from  Common 


Contents  of  the  labeled  commons  /NTPS/,  /XYZ/  and  /ZYZ/. 


*See  Part  I,  Vol  I 
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Data  from  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

NSTRIP 

NS1ZB(50,  2) 

JL1M(50,  3) 

XIS LCT(50,  7) 

j 

NTP10 

See  2.  5. 3 

CK(4,  4,  50) 

NTP9 

See  2.  5. 1 9 

Data  to  Common 

None 


Data  to  T apes 

Usually  none.  However,  in  special  cases,  when  dynamic  stability 
derivatives  are  also  required  (input  flag  NPR3=3),  tape  unit  2  is  used  for 
saving  the  total  lift-and  moment  coefficients  computed  for  the  steady  case 
(KR(1)  =  0.0).  Subsequently,  this  information  is  read  from  tape  unit  2  and 
is  used  in  the  calculation  of  the  dynamic  stability  derivatives  for  the 
unsteady  case  (KR(2)  >  0.  0)  -  see  2.  4.  2,  in  Part  I,  Vol  I, 


Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION 

BR 

br>  c/2 

\ 

32 

b^*  s 

\  See  2.  5.  19 

KR 

k 

r 

) 

ACAP, 

REFARA 

A,  S 

Reference  area  -  input 
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MNEMONICS 

SYMBOL 

DESCRIPTION 

NRF 

Number  of  reduced  frequencies 
for  case  -  input 

NSV 

Number  of  strips  on  vertical  panels 
(fin)  of  T-tail  cases;  0  otherwise 
-  input 

NYAW 

Yaw  flag  -  input 

IMX 

Number  of  rigid  body  modes 
for  AlC-type  case-from  GENB 

NPC 

Camber  flag  I 

NDELT 

6 

Symmetry  flag  )  -  input 

NPR3 

Print  flag  1 

2.5.21  SUBROUTINE  FINAL  (BR.  B2.  ACAP.  NSV,  NBV.  IMX. 

NYAW,  NPC.  NDELT,  FGAMMA) 


Functional  Description 

Subroutine  FINAL  is  used  for  AlC-type  cases  only.  It  reads 
the  solutions  for  the  AlC-type  right  hand  sides  from  tape  NTP3  and  computes 
the  pressures  for  the  rigid  body  modes  for  all  boxes  of  ail  lifting  surfaces. 
These  are  printed  along  with  the  box  number  and  the  fractional  chordwise 
location  of  the  load  points  (1 /4-chord  of  box).  Subroutine  FINAL  can  be 
bypassed  by  setting  the  data  flag  NPR3  to  0. 

Data  from  Cards 


None 


Data  from  Common 

Contents  of  the  labeled  commons  /NTPS/,  /XXZ/,  and  /YZY/ 
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Data  from  Tapes 


MNEMONICS  SYMBOL  TAPE 


NSTRIP 
NSIZB(50,  ?■) 
LIM(50,  3) 
CPM(400,  4) 


TAPE 

DESCRIPTION 

|  NTP10 

See  2.  5. 3 

NTP3 

Solution  matrix  for  the 

'  „ 

AlC-type  right  hand  sides 

Data  to  Common 


Data  to  Ta 


Data  in  Argument  List 


MNEMONICS 

h  ■■■ 

SYMBOL 

DESCRIPTION 

BR 

b 

r 

J 

B2 

%•  s 

!  See  2.5.  20 

ACAP, 

REFARA 

A,  S 

NSV 

' 

NBV 

Total  number  of  boxes  on  vertical 
panels  (fin  )  of  T-tail  cases  -  input 

IMX 

\ 

NY  AW 

(  See  2.1.2 

NPC 

\ 

NDELT 

6 

FGAMMA(50) 

Array  of  dihedral  angles  for  all 
strips  of  all  lifting  surfaces. 

2.  5.  22 


SUBROUTINE  GENF  (NB0X,  NT0T,  NMD.  NDELT.  NST. 

LIM.  FL.  XCAP,  XU,  EV.  SDELX. 

B2.  Y.ZZ.XOC.  NBV) 

Functional  Description 

Subroutine  GENF  is  used  for  modal  type  cases  only.  It  reads 
the  solutions  for  all  modes  from  tape  NTP3  and  writes  the  pressures  for 
all  modes  and  all  boxes  along  with  the  x,  y,  z  coordinates  of  the  boxes. 

It  also  computes  the  sectional  lift  and  moment  coefficients  for  all  strips 
of  all  lifting  surfaces  as  well  as  the  total  lift-  and  moment  coefficients  - 
see  Eqs.  (2.  4-6)-(2.  4-1  ^Subroutine  GENF  also  computes  the  generalized 
forces  for  all  pressure  modes  and  deflection  modes;  see  Eq.  2.4-5* 

Data  from  Cards 


None 


Data  from  Common 

Contents  of  the  labeled  common  blocks  /NTPS/,  /XYZ/  and  /B0DY/. 


Data  from  Tapes 


MNEMONICS 

SYMBOL 

TAPE 

DESCRIPTION 

P(400) 

N 

NTP3 

Solution  matrix  for  modal 
type  data  read  in  column 
order.  Dimension  of  matrix 
is  NB0XxNMD. 

P(400) 

I  AC  (f)j 
p  ! 

NTP7 

Matrix  of  pressures  for  all 
body  elements  and  modes  read 
in  column  order.  Dimension 
of  matrix  is  NBExNMD.  This 
data  is  read  for  body  cases  only 

B(400) 

NTP10 

Integration  matrix  for  modal 
type  data  read  in  column  order. 
Dimension  of  matrix  is  NT0TxNMD 
where  NT0T  =  NB0X  +  NBE. 

Data  to  Common 


None 

*See  Part  I,  Vol  I 
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Data  to  Tapes 


Data  in  Argument  List 


MNEMONICS 

SYMBOL 

DESCRIPTION 

NDELT 

NB 

6 

Symmetry  flag  \ 

Number  of  bodies  \  input 

NSTRIP 

Number  of  strips  ) 

NB0X 

Total  number  of  boxes  on  all  panels 

NT0T 

n 

NB0X  +  NBE,  where  NBE  =  total 
number  of  body  elements 

NBV 

\  T -tail  input 

NSV 

t 

NMD 

NM,  m 

N'.  mber  of  modes  -  input 

LIM(50,  3) 

See  Sec.  2.  5. 3 

ACAP, 

A.  S 

\ 

RE FAR A 

\ 

FL, 

c,  2b 

1 

REFCHD 

r 

1 

B2, 

s,  b/2 

f 

REFSPN 

} 

EV  (400) 

e 

1  See  Blank  Common,  Sec.  2.  3.  1 

Y(400) 

y 

I 

ZZ  (400) 

Z 

J 

SDELX(400) 

Ax 

s 

XIJ  (50) 

Xle 

\  See  Common  /XXZ/,  Sec.  2.  3.2 

X0C(4OO) 

£/c 

J 

2.6  BLOCK  DIAGRAM 


The  basic  block  diagram  for  the  computer  progrm  H7WC  is  shown 
on  the  next  four  pages.  All  logical  steps  in  the  "main"  part  of  the  program 
are  indicated,  and  the  basic  functions  of  all  subroutines  that  are  called 
from  "main"  are  outlined.  A  complete  description  of  all  subroutines  — 
including  those  that  are  not  called  directly  from  "main,"  and  therefore 
are  not  indicated  in  the  block  diagram  —  can  be  found  in  Sections  2.  5.  1 
through  2.  5.  22. 


compute  rigid  body 
mode  pressures 


2.  7  PROGRAM  LISTING 


The  following  pages  (ppv.  95  -  143  )  contain  the  complete  listing  of 
program  H7WC.  The  individual  subroutines  of  the  program  (including  the 
"main"  part  of  it)  are  listed  in  the  order  of  their  description  in  Section 
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